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A QUANTITATIVE THEORY OF GENETIC RECOMBINATION 
AND CHIASMA FORMATION® 


R. A. FiIsHER 


Cambridge University 


1. INCOMPLETENESS OF GENETICAL THEORY 


I have undertaken to talk on the quantitative theory of the forma- 
tion of chiasmata and of genetic recombination with some diffidence, 
because the ideas I shall try to put before you are so new as to be in 
some respects immature, or half developed. They have arisen in con- 
nection with the experimental work with mice carried out in the Depart- 
ment of Genetics at Cambridge during the last few years. I can reason- 
ably attempt in the course of one lecture to give the background of 
reasoning on which these ideas are based, but there is one difficulty 
which I cannot hope to overcome, and that is to give the detailed mathe- 
matical developments, largely due to my colleague George Owen, by 
which the primary scheme is connected with its detailed and experi- 
mentally verifiable consequences. 

The subject of crossing over has for many years seemed to me to be 
that part of genetics in which teachers are furthest from having a com- 
plete theory to present. Let me illustrate this as simply asI can. We 
may think of a strand or chromosome in a gamete, having on it two 
points or loci marked by genetic factors. Then, in the preceding 
meiosis the number of breaks in the intervening tract, involving inter- 
change of parts of the parental chromosomes, is one or another of the 
numbers 

These events are mutually exclusive, and if their respective 


P * Lecture delivered September 5, 1947, before the First International Biometric Con- 
erence. 
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probabilities are represented by : 


Poy Pr» Pa, Day 5 


then, since one or another of these events must have occurred, it follows 
that the sum of their probabilities is unity, or that 


Pot Pi, +Pet+Pst...-.=1. 


That is simple and obvious, but I need not remind geneticists that 
this simple specification of the rules of gamete formation is not the 
usual way in which the subject has been discussed. It requires no 
mention to be made of the possible configurations of the four-strand 
stage. It specifies, however, all that is relevant to the purely genetic 
question :—what kinds of offspring will be produced by any given kind 
of parent. The genetically observable, or quasi-observable, quantities 
can in fact be specified in terms of the series of p values. 

First, the recombination fraction is directly observable by breed- 
ing tests; it is evidently equal to the sum of the probabilities of an 
odd number of breaks. I.e.— 


Next, the map distance, z, may be defined as the expected number 
of breaks in the tract, so that 


r= 


These relations may be expressed somewhat more clearly by defin- 
ing a second series of quantities s, to represent the frequencies of 
r or more interchanges in the segment, so that 


= Pi, Dr = — 


Then, the equations, 
Y 


express the two genetic observables in terms of the s-series. 

_ Now, any adequate theory of crossing over must supply the means 
of expressing the entire series of p-values, or of s-values, in terms of the 
positions of the genetic markers on the chromosome, and their relation- 
ships to its centromere and to its ends. How far genetical theory is in 
fact from meeting this requirement may be illustrated by the attempts, 
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of obvious practical value, to express z in terms of y. Since only y is 
directly observable this is a primary need in making a linkage map. 
Early attempts at this problem were made by Haldane, and a few years 
ago a remarkably effectual formula was put forward by the Indian 
mathematician, Kosambi, namely that 


2y = tanh (2z). 


Kosambi’s formula leads to an addition theorem for the recombina- 
tion fraction of the sum of 2 segments, 


Y2 
1+4y.y2’ 
and with three genetic markers this may be tested by direct comparison 
with the data. In such tests Kosambi’s formula usually shows up very 
well, yet it can scarcely be accurate in all cases, for it takes no account 
of the position of the segment in relation to the centromere, a factor 
which has long been known greatly to influence interference. 

In the second place, even if it were accurate in expressing a fixed 
relationship between z and y, a formula of this kind could not supply 
a complete theory. For work with three factors the relationship be- 
tween z and y is sufficient, but with four or more factors segregating 
simultaneously it is easily seen to be insufficient. With four factors, 
for example, there are eight pairs of complementary genotypes, and 
seven constants are needed to specify their relative frequencies; yet an 
addition formula, such as that given above, can only give the recom- 
bination fractions between the six pairs of loci. A further parameter 
is needed to specify the gametic series. In fact, only the complete 
series of expressions for s:, $2, 8; . . . in terms of some common pa- 
rameters can give a theoretically complete aceount of recombination. 


= 


2. MATHER’S SKETCH OF CROSSING OVER AS A RANDOM PROCESS 


I should like to recall at this point some quite early proposals made 
by K. Mather in which the formation of chiasmata in a four strand 
pair of chromosomes was regarded as a random process, starting from 
the centromere, in such a way that the interval between the centromere 
and the chiasma next to it was assigned a definite frequency distri- 
bution, and the interval between the first chiasma and the second was 
treated in the same way. It is evident that with proper specification 
of these distributions, calculation should be able to provide the fre- 
quencies of all definable configurations, specifying the incidence of 
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simultaneous sets of chiasmata in the same chromosome: 

I have followed Mather’s ideas in this matter with some variations. 
First, I consider only the configurations provided by sets of chiasmata 
on the same strand, for these are enough to specify the genetic situa- 
tion. Next, it appears that good agreement with the facts may be ob- 
tained by introducing a metric which is not to be identified with map 
distance, or with physical length, in terms of which interference is a 
constant property at all parts of the strand. In this metric it appears 
to be sufficient, at least as a good approximation, to treat the cen- 
tromere as if it were an obligatory chiasma. In terms of this metric 
it is then theoretically possible to caleulate the series s,, s2, s3,.... 
for any chosen points, and in particular the values of x and y, so as to 
give full information about recombination for all regions of a given 
map. 


3. THE SEX CHROMOSOME OF THE HOUSE MOUSE 


Some years ago, when Mather was working with me at the Galton 
laboratory, we carried out, with the able assistance of Mrs. S. Holt, an 
experiment which at the time must have seemed rather extravagant. 
Although the laboratory rodents had long been the best studied species 
among the mammals, yet sex linkage had not appeared among them. 
In these circumstances, it seemed to us worth trying whether any of the 
factors we had readily available were, as was possible, carried in the 
pairing segment of the sex chromosome, and in consequence exhibited 
partial sex linkage. In the event, none of the six factors we tested: 
pied, agouti, brown, dilute, wavy (wv,), and light head, appeared to be 
linked with sex in a test with about twelve hundred mice. Yet before 
our experiment was ready for publication, J. B. S. Haldane, working 
in the same College had made a search for partial sex linkage in Man 
and had found pedigrees of several rare defects appearing to be so 
inherited. The test tor a similar phenomenon in mice, was evidently 
not so hopeless. 

When in 1943 I moved to Cambridge I was able, through the gen- 
erosity of The Rockefeller Foundation and the kindness of Dr. Snell of 
the Roscoe B. Jackson laboratories to increase the variety of my mouse 
stocks, and, about two years ago Mrs. M. Wallace, then Miss Margaret 
Wright, noticed a certain irregularity in the sex ratios of some of the 
lines which we kept segregating for the linked factors wavy-two and 
shaker-two. In consequence we decided to set up an experiment which 
we hoped would be competent to discover what was really happening 
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at these two loci. The effect to be demonstrated was evidently not very 
large, and there were reasons to fear, as there usually are with reces- 
sive factors, and especially shakers, that the evidence would be dis- 
turbed by differential viability. Consequently, the experiment had to Ses 
be carefully designed. ae 
There are four possible types of male heterozygous in these 2 fac- i 
tors, for the animal may have received from his father, either wavy or 
shaker or both or neither, with a complementary set of genes from the 
mother. Of these four male types that in which neither recessive 
comes from the father and both from the mother, i.e. 


‘wv, sh, X/++Y (Tricoupling) 


may be regarded as in coupling for all three factor pairs, while the 
other three types differ from it respectively by interchanging sex (X, 
Y), wavy (+, wv,), and shaker (+, sh.) giving therefore 


++X/wv, sh, Y 
+ sh, X/wv.+Y 
wv, + X/+ sh, Y 
Each of these male types mated with homozygous wavy shaker does, 
wv, sh, X/wv, sh, X 


produces eight genotypes, four male and four female, which fall into 
four pairs of complementary genotypes 


Wavy shaker females and normal males A 
Normal females and wavy shaker males B 
Shaker females and wavy males Cc 
Wavy females and shaker males D 


The pairs of genotypes B, C, D are derived from the pair A, by the 
same three interchanges of male with female, of normal with wavy, and 
of normal with shaker. 

In relation to any one of the types of heterozygous fathers these 
four pairs of genotypes are produced by four modes of gamete forma- 
tion of which one, in the case of three loci in the same chromosome, may 
be identified as old combinations, while in the others recombination has is 
taken place in one segment or the other, or in both. Without knowing pa 
the order of the loci we cannot tell which of these is which; but in all. 
cases the three modes of formation involving recombination will differ 
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from the first, in which no recombination has occurred, by the three 
operative interchanges for sex, wavy and shaker. 

Consequently the data obtained when all four male types are used 
will consist of sixteen frequencies, each of a pair of complementary 
genotypes arranged in a symmetrical four by four Latin square, as 
shown below: 


Type of gamete formation 


No Sex Wavy Shaker 

Type of male change change change change 
wv, sh, X/++Y A B Cc D 
++X/wv, sh, ¥ B A D “4 
+ 8h, X/wv, + Y C D A B 
wv, + X/+ sh, Y D Cc B A 


The experimental purpose of obtaining data in this form is to allow 
of the estimation of the frequencies of the four modes of gamete forma- 
tion, free from disturbances introduced by unknown differences in the 
viabilities of the four pairs of genotypes, and of arbitrary differences 
in the numbers of offspring bred to the four types of male. This sta- 
tistical problem is intricate and has never been adequately discussed. 
Fortunately, however, the particular data obtained with the factors 
wavy-two and shaker-two in mice offer no difficulty, for there were in 
our results no signs of differential viability, as appeared from the com- 
plete homogeneity of the four distributions shown by the rows of the 
table. 

The relative frequencies of the four modes of gamete formation may 
in such a case be estimated directly from the total frequencies in the 
columns. Out of about 450 mice bred it was found that the first 
column contained about 29%, the second 41%, and the remaining 
columns about 15% each. I need not emphasize that such a result ap- 
peared highly paradoxical. For, in the first place, if sex linkage were 
not involved, we should expect equality in the first and second columns, 
as we would also in the third and fourth ; whereas the first two columns 
differ consistently for all types of male, and in the aggregate by over 
50 mice. Next, if sex-linkage were involved we should expect entries 
either in the third or in the fourth column to be very rare double cross- 
overs, but these two classes are nearly equally frequent. Finally with 
incomplete sex-linkage we should expect the sex interchanges to be less 
numerous than the cases of no recombination. 
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In respect to recombination frequencies we have :— 


More accurately 
Wavy-shaker 15+15=30% 31.1% +2.17 
Sex-shaker 41+15=56% 56.1% + 2.33 
Sex-wavy 41+15=56% 56.7% + 2.33 


So that the paradoxical nature of our data may be restated by saying 
that both the locus for wavy, and that for shaker, have significantly 
more than 50% recombination with sex ; moreover, these recombination 


fractions are nearly equal, although the loci must be 30 to 35 map units 
of distance apart. 


4. TOPOGRAPHICAL INTERPRETATION OF THE DATA 


In the male mouse it is to be supposed that the inheritance of sex 
takes place as if a sex determining locus were situated at the point of 
separation of the pairing segment from the differential segment. For 

Figure 1 


Frequency Distribution of Length of 
Intercept Between Successive Breaks 


RELATIVE FREQUENCY 


0.5 1.0 15 2.0 
METRICAL LENGTH OF INTEREST, 


purposes of calculation it is convenient to take this point as origin. 
The occurrence over a considerable length of the pairing segment of 
recombinations with sex exceeding 50 percent, becomes more intelligible 
if we consider that the complete absence of chiasmata in the differential 
segment of these strands must mean an absence of interference from 
this side; so that, on a metric chosen for uniformity in respect of inter- 
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ference, this region will be favored by chiasmata. Moreover, calcula- 
tion taking any plausible form (see fig. 1) for the distribution 
of the intercept length such as assigning to it the frequency element 


sech tanh ) 


where wu is the metrical value, and the coefficient ; 7 is introduced to 


make the metric conform to map distance in regions far from the arm 
ends, or from the centromere, shows that this expected effect is realized. 

If now numerical values of s,, S2, $3, .. . . are tabulated for various 
distances from the end, the value of the sex recombination fraction, y, 
and the map distance, z, may each be obtained in terms of wu. Then, 
plotting the recombination fraction y against the map distance zg, it 
appears that y rises at first in equality with z, but passes the 50 percent 
value at a map distance of about 60 map units and rises to a long flat 
maximum between 55% and 56% recombination. 


Figure 2 
Recombination Fraction in Terms of Map Distance 
% 
50a 
= 
< 
Zz X- 
3 
104 


MAP, DISTANCE 


The genetical findings with wavy and shaker would therefore be | 
predicted by theory if these two loci were situated about 30 units apart 
in the neighborhood of the maximum recombination fraction and 
probably on either side of it. 

In a more detailed study George Owen has shown how the effects of 
finite lengths of the chromosome arms may be taken into account; from 
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his work it appears also that recombination fractions exceeding 50 per- 
cent are characteristic of a wide class of intercept distribution func- 
tions, provided interference is sufficiently pronounced. He has found 
models in which explicit analytic functions may be used to replace the 
somewhat cumbrous numerical integration I had employed. 

As I mentioned in my opening remarks, very much detailed work 
remains to be done, although, so far as we have explored the question, 
the genetical data have shown themselves entirely conformable with the 
numerical details of the theory. The exact degree of confirmation ob- 
tainable from the available data for Drosophila will, however, not be 
clear until a more thorough study has been made of the relationship 
of coincidence values, appropriately estimated, to the chromosome 
regions in which they are found. 


DISCUSSION 


D. G. Catcheside. The occurrence of recombination in excess of 
50 per cent requires that there should be chromatid interference favor- 
ing a frequency of four-strand doubles in excess of that expected in the 
absence of chromatid interference. 

When two successive chiasmata occur in an arm of a pair of chro- 
mosomes, the relationships between the four chromatids may be such 
that two strands (complementary), or three strands (diagonal) or four 
strands (reciprocal) take part in the crossing-over at the two chiasmata. 
If the participation of a given strand in one chiasma does not affect the 
chance of its participation in the other chiasma in the arm, the strand 
relationships should occur in double crossing-over with the frequencies 
of a quarter complementary, a half diagonal and a quarter reciprocal. 
This absence of chromatid interference is found in Drosophila melano- 
gaster and with such absence the amount of recombination between 
linked genes cannot rise about 50 per cent. 

Chromatid interference may be of two types; either (a) the par- 
ticipation of a chromatid in one chiasma reduces its chance of partici- 
pation in a second chiasma (positive interference) or (b) participation 
of a chromatid in one chiasma increases its chance of participation in 
a second chiasma (negative interference). The latter, in which the 
upper limit of recombination would be less than 50 per cent, has been 
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observed in Neurospora crassa, where 2-strand doubles are in excess 
and 4-strand doubles are deficient. The-former type, which could in 
some circumstances lead to an amount of recombination between two 
loci in excess of 50 per cent, has not so far been observed. It may be 
that crossing-over in the heteromorphic pair of sex chromosomes in the 
male mouse is of this type, but it must not be assumed that positive 
chromatid interference is of general occurrence. 

While it is true that most genetical experiments have shown an 
absence of chiasma interference across the centromere, recent work of 
Callan and Montalenti (Journal of Genetics 48, 119; 1947) has con- 
firmed the occurrence in Culex pipiens of positive chiasma interfer- 
ence across the centromere as previously recorded by Patan. Thus the 
centromere cannot always be taken as a null point. 


Hilda Geiringer. Professor Fisher’s new and extremely interest- 
ing theory of recombination suggests the following question: The com- 
plete linkage of m factors is characterized by M = 2"*-1 parameters. 
In fact, denote the constitution of an organism by 22, . . . 2m; Y1, 
Yo, ... Ym) Where the maternal (paternal) heritage with respect to 
factor is denoted by x; or y; respectively. A gamete may then be 
described, for example, by 21, Z2, Ys, Ys; Ys, - - - Xm; there are 2" gametes 
of this sort and there exists a probability distribution corresponding to 
these 2™ possibilities. Admitting obvious symmetry relations this num- 
ber, 2”, will be divided by two and we thus have 2" probability values 
with sum equal to one, hence M parameters. This ‘‘linkage distribu- 
tion’’ (1.d.) is completely equivalent.to the cross-over distribution, or 
‘*recombination distribution’’ (r.d.) which (in the ‘‘linear theory’’) 
specifies the probabilities that in each of the (m—1) intervals between 
the loci of linked factors recombination takes place or does not take 
place. 

On the other hand this r.d. is, in general, not equivalent to the 
N=m(m-—1)/2 recombination values. In fact M>WN for m= 4. 
From tlie r.d. (or from the |.d.) the recombination values follow, but 
not vice versa. It is, however, an aim to be able (by means of observa- 
tion and theory) to reduce the number M of independent parameters 
and to express the M values by means of the N recombination values 
or possibly by part of them, e.g. by the (m—1) recombination values 
between adjacent factors. May I ask Professor Fisher’s opinion regard- 
ing this question, in the light of his new theory? 
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Alexander Weinstein. Recombination exceeding fifty per cent. 
The recombination frequency between two loci depends on (1) the 
number of chromatids that cross over at a level, (2) which chromatids 
these are (whether homologous or sister strands), (3) the number of 
levels at which crossing over occurs, and (4) the relation of exchanges 
at different levels to one another (whether regressive, progressive, or 
digressive). 

If at any level only two chromatids of a tetrad cross over and these 
are always homologous and not sister strands, then, on certain assump- 
tions, in crossover tetrads with n exchanges between two loci, the pro- 
portion of chromatids showing recombination between these loci is 4 
when n is odd and 4 — 4 (F — H)”? when n is even; where F is the chance 
that an exchange is regressive and H the chance that it is digressive 
with respect to an adjacent exchange. These formulas apply not to 
chromatids from all tetrads taken together, but to chromatids from 
tetrads of a specified rank n. 

If F = H, then in a group of crossover tetrads of rank n, the recom- 
bination frequency of emerging chromatids remains 4 as nm increases. 
If F exceeds H, the frequency will oscillate between 4 and a lesser 
value which approaches $ as n increases indefinitely. If F is less than 
H, the recombination frequency oscillates between values greater and 
values less than 4, and these values approach 4 as 7 increases indefi- 
nitely. That is, the recombination frequency in chromatids derived 
from tetrads of rank n exceeds 4 if F is less than H and 5 is odd ; and 
if such tetrads are sufficiently numerous, the recombination frequency 
in chromatids from all tetrads taken together may exceed 50 per cent. 

These results assume that F and H are constant throughout the 
chromosome. If either or both vary according to the regions involved, 
the recombination frequency in chromatids derived from tetrads with 
exchanges in n specified regions between two loci remains 4 when n is 
odd; and when n is even, the recombination frequency becomes 
... where the subscript 
indicates the second of the two adjacent regions involved ; for example, 
F, is the chance that an exchange in region 4 is regressive with respect 
to an exchange in region 3. The recombination frequency in chro- 
matids derived from a crossover tetrad of rank n would exceed 4 if F 


1 Weinstein, Alexander, “Mathematical study of multiple-strand crossing over and 


coincidence in the Sacrnecomes of Drosophila,” American Philosophical Society Year 
Book 1937, pp. 227-228, 1 
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is less than H in an odd number of factors, provided that in the other 
factors F exceeds H. If F equals H in one or more factors, the recom- 
bination frequency is 3. 

The recombination frequency might also exceed 50 per cent if an 
exchange between two chromatids of a tetrad were accompanied by an 
exchange between the two others at the same level. This would yield 100 
per cent recombination in chromatids derived from tetrads with cross- 
ing over at one level: the result would be the same as if the two ex- 
changes were at different levels and digressive with respect to each 
other. If two exchanges at the same level are treated as mutually 
digressive, the formulas already given can be applied, whatever the 
number of levels at which there is crossing over; provided that n is 
taken as the number of exchanges, which is now no longer identical 
with the number of levels at which exchanges occur. 

The relative frequency of regressive, progressive, and digressive 
crossing over also affects other genetic functions, including coincidence 
of various types. In attached X-chromosomes, the frequency of homo- 
zygosis in individuals derived from tetrads of rank 2 or higher depends 
on the frequency of progressive exchanges, G, which equals 1-— F — H. 
Thus in suitable material the values of F and H calculated from the 
recombination frequencies can be subjected to additional tests. 

(Such a notation was presented orally, but is omitted from the 
written proceedings subject to publication elsewhere.) 

J. Lederberg. Professor Fisher has given us a strong justification 
of what he terms a genetic as opposed to a cytogenetic approach to the 
problem of interference in crossing-over. We shall, however, all be 
interested to see in the further development of this analysis the extent 
to which his genetic formulation is consistent with, or even relevant to, 
cytogenetic data on interference in the four-stranded bivalent in which 
the physical phenomenon of crossing over is believed to occur. 

In bridging this gap, the usual procedure is the construction of 
many tetrad diagrams in which the effects on the gametic output of 
combinations of crossovers involving various strands at different levels 
are worked out geometrically. The geometry is, however, sufficiently 
simple to be amenable to convenient notation, by the use of which the 
effects of crossing over can be elucidated algebraically. An algebra is 
most useful perhaps in respect to the formulation of theories in terms 
of restraints on the values of its terms, which may be more easily 
manipulated than restraints geometrically expressed. 
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R.A. Fisher. In reply to Mr. Catcheside, I am surprised and some- 
what taken aback that he should be of the opinion that chromatid inter- 
ference is absent in Drosophila melanogaster. Perhaps he only means 
that the interference observed could conceivably be explained without 
specifically chromatid interference ; but this seems to be quite another 
matter. It is at least obvious that the Drosophila data accords excel- 
lently with the theory here under discussion. 


In reply to Dr. cemrincer: The new theory does supply formulae, 
in some cases evaluable easily, and in others with more difficulty, 
which give the frequencies of all possible recombinations to be 
found in gametes from a multiple heterozygote. These are given in 
terms of the metric-values of the different segregating loci relative to 
the centromere, of the arm lengths of the two parts of the chromosome, 
and, possibly, of an additional parameter available for — the 
intensity of interference. 


In respect to the points raised by Dr. WEINSTEIN and others in 
the discussion, I should have made it clear that my theory does not 
depend on the whole complexity of the configurations possible in the 
four-strand stage, but only on the frequency distribution of the dif.- 
ferent configurations possible to a single strand received by a gamete. 
This simpler distribution is sufficient for the purely genetic purpose of 
predicting the frequencies of the different genotypes derived from a 
multiple heterozygous parent. 
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LA RELATION D’ALLOMETRIE 
SA SIGNIFICATION STATISTIQUE ET BIOLOGIQUE* 


GrorGES TEISSIER 
Station Biologique de Roscoff 


Le chapitre de 1’étude de la croissance dont les progrés récents ont 
été les plus remarquables est sans contredit celui qui traite des change- 
ments de proportions des différents organes ou, autrement dit, de la 
croissance relative. Les recherches sur cette question, fort nombreuses 
depuis une vingtaine d’années, ont fourni matiére a plusieurs mises au 
point, inégalement développées.* 

La question est actuellement si vaste et si complexe, qu’il apparait 
qu’en un exposé limité, il faut se contenter, ou de ne traiter que les 
grandes lignes du sujet, en laissant de cété, non seulement bien des 
points importants, mais aussi la plupart des difficultés qui y figurent 
encore, ou de s’attacher 4 une question beaucoup plus limitée, en essay- 
ant de l’approfondir. C’est 4 ce dernier parti que je me suis arrété 
et il ne sera question ci-aprés que d’un probléme unique, qui, 4 la 
vérité, se trouve au centre méme du sujet. Dans ce qui va suivre, nous 
allons chercher la signification exacte de la relation d’allométrie qui, 
depuis plus de vingt ans est |’instrument indispensable de toute étude 
sur la croissance relative. 

Le fait fondamental en la matiére est 1’existence d’un type de rela- 
tion qui, convenablement ajusté, permet de décrire exactement les 
croissances relatives les plus diverses. Chez tous les animaux appar- 
tenant 4 tous les grands embranchements étudiés jusqu’a présent, les 
relations entre, d’une part, les dimensions ou le poids d’organes trés 
différents par leur nature, leur structure ou leur fonctionnement et, 
d’autre part, les dimensions ou le poids du corps tout entier, ou d’un 
organe pris comme référence, peuvent étre traduites par des équations 
de la forme: 

y = bs 


‘ * Lecture delivered September, 1947, before the First International Biometric Con- 
‘erence. 


1G. Teissier, Travaug Station Biologique de Roscoff, 1931, 9; J. S. Huxley, Problema 


of Relative Growth, London, 1932; G. Teissier, Act. Sc. Hermann 95, 1934; Annales 


Physiologie et Physico-Chimie Biologique, 12, 1936, et Act. Sc. Hermann, 1936; E. C. R. 
Reeve et J. S. Huxley; P. B. Medawar; O. W. Richards et A J. Kavanagh, Essays on 
Growth and Form, Oxford, 1945. 
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Fig. 1 
Allométrie de croissance du céphalothorax de Carcinus moenas. Le croquis 
donne la définition des dimensions dont les mesures sont reportées sur le graphique. 
Coordonnées logarithmiques. g=1,00 pour a a; q=0,98 pour b b; q=0,93 pour 
ec; a=0,86 pourdd;a= 0,74 pour e e. 
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ou } et a sont des constantes. En coordonnées logarithmiques, ces ares 
de courbe puissance se traduisent par des droites de pente a: 
log y = log b+ a log z 

Lorsque la croissance d’un organe est plus rapide ou moins rapide que 
celle du reste du corps, ou de l’organe de référence, et que la relation 
précédente est vérifiée, nous disons qu’il y a allométrie. Sia > 1, 
l’allométrie est dite positive ou majorante ; si a < 1 elle est négative ou 
minorante. Dans le cas od a= 1, la croissance de l’organe et celle du 
corps sont paralléles: il y a tsométrie. 

Des deux paramétres qui définissent la relation d’allométrie, le plus 
important a, pente de la droite figurative du phénoméne en coordonnées 
logarithmiques, sera la constante actuelle d’équilibre; b qui représente 
la taille de l’organe y lorsque z=1 a la méme signification que les 
indices dont font usage les anthropologistes : ce sera l’indice origine. 


Avant de commencer un examen approfondi de cette relation, il n’est 
pas inutile, pour justifier le temps que nous lui consacrerons, de rap- 
peler l’importance du réle qu’elle joue actuellement dans 1’étude de la 
croissance relative et dans la biométrie générale. 

En premier lieu, nous devons insister sur le fait qu’elle traduit en 
général avec la plus grande fidélité les modalités de la croissance 
moyenne. On pourrait le prouver par plusieurs centaines d’exemples 
également probants. L’étude statistique des cas pour lesquels on 
posséde des mesures suffisamment nombreuses permet par ailleurs de 
démontrer que les écarts qui se produisent de part et d’autre de cette 
loi moyenne peuvent étre considérés comme fortuits. Elle permet en 
outre de comparer, avec autant de sécurité, des animaux adultes appar- 
tenant soit 4 la méme espéce, soit 4 des espéces différentes. C’est a 
ee titre qu’elle est devenue un instrument trés efficace dans 1’étude des 
séries phylogéniques ou dans la recherche des affinités des races d’une 
méme espéce ou des espéces d’un méme genre. 

Certains organes ou certaines parties du corps, dont la croissance 
ne peut étre représentée par une seule courbe puissance, suivent suc- 
cessivement plusieurs relations d’allométrie. Sur les graphiques log- 
arithmiques les droites peuvent se raccorder par un point anguleux, ou 
étre séparées par une discontinuité plus ou moins ample, qui correspond 
& un stade oti le poids (ou la longueur) de l’organe n’est plus fonction 
du poids (ou de la longueur) du corps. Une courbe de croissance rela- 
tive peut ainsi présenter deux sortes d’accidents et se décompose dans 
les deux cas en troncons réguliers exactement définis. 
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Les points anguleux et les discontinuités qui jalonnent les courbes 
de croissance des différents organes d’un méme animal ne sont pas 
répartis de fagon queleonque, mais apparaisscnt au contraire simul- 
tanément dans des organes divers. La croissance d’un animal se divise 
ainsi naturellement en étapes distinetes. Pendant chaecune de ces 
étapes, les divers organes et les divers constituants biochimiques erois- 
sent, l|’un par rapport 4 l’autre, en suivant des relations d’allométrie 
simples, la composition chimique, la structure et la forme variant de 
fagon continue. Les étapes successives sont séparées l’une de 1|’autre 
par des stades critiques de durée assez bréve ou, sans que le poids ou la 
taille augmente sensiblement, les régles—qui présidaient jusque 1a a 
la croissance changent brusquement. Ces stades critiques se marquent 
sur les diverses courbes de croissance relative par des points anguleux 
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Fig. 2 
Etapes de la croissance du propodite de la pince chez les Maja squinado miles. 
Longueur du propodite en fonction de la longueur du Céphalothorax. Coord mnées 
logarithmiques. Male: q=1,09 pour la premiére étape; a=1,35 pour la deus iéme; 
a= 1,93 pour la troisitme. Femelle: a= 1,09 pour la premiére étape; a=1,1¢ pour 
la deuxiéme; a=1,16 pour la troisiéme. 
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ou des discontinuités. Du point de vue biochimique ils correspondent 
a des changements plus ou moins profonds dans 1’équilibre réciproque 
des divers constituants de l’organisme. Du point de vue physiologique, 
ils correspondent fréquemment 4 des modifications importantes dans 
le jeu des corrélations humorales. 

Cette conception de la croissance relative, que j’ai développée autre- 
fois me parait, encore aujourd’hui, correspondre assez exactement a 
la réalité bien que sur certains points, notamment en ce qui concerne le 
role des stades critiques, l’accord des spécialistes soit loin d’étre fait. 
Il n’est pas question de la défendre ici mais, qu’elle soit exacte ou 
qu’elle soit fausse, son succés définitif, ou son échee, dépendront en 
définitive de la validité plus ou moins grande de la relation d’allométrie. 
Aucune des interprétations physiologiques, voisines d’ailleurs 1’une de 
l’autre qui ont pu étre suggérées par Huxley, par Robb, ou par moi- 
méme, n’a jusqu’ici recueilli l’assentiment général. Peut-étre en 
trouvera-ton une plus satisfaisante, mais il restera encore 4 expliquer 
que la méme formule convienne 4 tant de cas si disparates et 4 com- 
prendre sa signification générale. 

Avant d’aborder ce probléme, une remarque préalable est indispens- 
able. La raison essentielle du succés de la notion d’allométrie est 
qu’elle ne fait pas intervenir l’age des individus dont elle compare les 
caractéristiques. C’est grace a cette circonstance qu’ont pu étre faites 
des études trés précises sur la croissance d’animaux qu’il est difficile, © 
ou méme impossible, d’élever au laboratoire. 

Une telle caractéristique n’est évidemment pas liée 4 la forme méme 
de la relation d’allométrie et toute autre relation ¢ (z, y) = 0 entre deux 
grandeurs mesurées sur un méme animal ‘‘courtcircuiterait le temps”’ 
aussi bien qu’elle. Elle postule, en revanche, la réalisation préalable 
d’une condition rarement énoncée d’une facon explicite, parce que 
sans doute trop évidente, mais qui est nécessaire 4 la réussite d’une 
étude sur la croissance relative. I] faut que, chez l’organisme étudié, 
la forme, la structure ou la composition chimique d’un individu placé 
dans des conditions données, ne dépendent que de ses dimensions et non 
pas de son age. 

Qu’une telle condition puisse étre réalisée au cours de la croissance 
n’est nullement une nécessité logique. Elle apparait comme une régle 
pratique qui se vérifie, rigoureusement pour de nombreuses espéces, 
approximativement pour beaucoup d’autres, pendant la plus grande 
partie de leur existence et qui, chez elles, ne se trouve éventuellement 
en défaut qu’a certains stades critiques du développement, ceci n’ex- 
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cluant naturellement pas la possibilité, pour quelques espéces, de 
répondre par un changement de forme 4 certaines modifications de leur 
milieu. Il est encore beaucoup trop tot pour tenter de dresser la liste 
des espéces appartenant a cette catégorie d’étres 4 forme stable, mais 
on doit cependant citer au premier rang d’entre eux le grand groupe 
des Arthropodes. Comme par ailleurs, ces animaux se prétent a des 
mesures multiples et généralement beaucoup plus précises que celles qui 
peuvent étre pratiquées sur les représentants des autres embranche- 
ments, ils constituent, avee certains Vertébrés, un matériel de choix 
pour l’étude des lois de la croissance relative. Il semble, en revanche, 
a en juger du moins par les recherches des auteurs japonais, que les 
Lamellibranches se prétent mal a ce genre d’investigation, leur forme 
dépendant assez largement a taille égale, non seulement des conditions 
de milieu, mais peut-étre aussi directement de 1’age. 

Nous laisserons de coté les animaux appartenant 4 cette deuxiéme 
eatégorie pour nous occuper uniquement de ceux chez qui, pendant une 
partie au moins de la croissance, deux lots d’animaux de méme taille 
et d’age différent ne sont pas discernables 1’un de |’autre et chez qui, 
réserve faite de la variabilité spécifique, la forme est fonction de la 
taille. 

Il existe une autre circonstance, toute différente en son principe de 
la précédente, ot il est légitime de comparer les dimensions x et y de 
deux organes en faisant abstraction du temps, ¢’est celle ot la com- 
paraison porte sur des animaux ayant cessé de croitre et chez qui, pour 
un temps, ou définitivement, forme et taille sont comme figées. Il en 
est ainsi pour les Mammiféres adultes et pour les imagos d’Insectes et 
1’on sait que c’est sur des animaux appartenant 4 ces groupes et rem- 
plissant ces conditions qu’ont été faites quelques-unes des premiéres ap- 
plications des relations d’allométrie. Celles qui portent sur les Mam- 
miféres essaient en général de décrire les relations d’allométrie pouvant 
exister entre les structures d’espéces voisines mais trés différentes de 
taille, tandis que celles qui portent sur les Insectes tentent d’analyser 
les différences de forme corrélatives des différences de dimensions a 
l’intérieur d’une méme espéce. Depuis, allométrie de taille et allo- 
métrie de croissance ont été étudiées, 4 l’aide des méme méthodes, par 
les mémes auteurs sans qu’aucun d’eux ait réussi 4 justifier entiérement 
V’emploi qu’il faisait d’une méme formule pour traduire des relations 
aussi dissemblables que celles qui peuvent unir un jeune animal et 
l’adulte correspondant, deux adultes de la méme espéce, ou deux adultes 
d’espéce différente. Encore moins a-t-on réussi 4 expliquer que la 
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relation d’allométrie puisse étre un instrument puissant d’analyse dans 
l’étude des variations intraspécifiques et interspécifiques de la com- 
position chimique des organismes. 

C’est 4 la recherche d’une interprétation générale de la relation 
d’allométrie que va étre consacré le présent exposé. 


I 


Plusieurs schémas physiologiques, assez voisins d’ailleurs l’un de 
l’autre, ont été proposés pour interpréter l’allométrie de croissance. 
Si aucun d’eux n’a, jusqu’a présent, recueilli ]’assentiment général des 
biologistes, il est permis d’espérer qu’un schéma analogue, mais plus 
parfait, sera reconnu quelque jour comme satisfaisant. I] est certain, 
en revanche, qu’un schéma mécaniste de ce type ne pourra s’appliquer 
qu’aux cas ot la comparaison porte sur les stades successifs du dével- 
oppement d’un méme individu et 4 ceux ot les individus comparés 
peuvent légitimement étre considérés comme fournissant l’image des 
différentes étapes d’une méme évolution. Il perdrait toute significa- 
tion si l’on voulait 1’appliquer aux cas ot la comparaison porte, en 
principe comme en fait, sur des individus différents et ot ia relation 
d’allométrie condense en une formule une infinité de possibilités dont 
chaque individu ne peut, dans les circonstances les plus favorables, 
réaliser qu’une seule. J’ai proposé en 1937, une interprétation statis- 
tique de cette allométrie de taille et j’ai pensé 4 cette époque qu’ il serait 
déraisonnable d’essayer de la transposer au cas de |’allométrie de crois- 
sance. Je crois cependant aujourd’hui a la possibilité d’une interpré- 
tation des phénoménes d’allométrie qui engloberait non seulement tous 
les faits morphologiques, mais aussi les faits biochimiques pour lesquels 
il n’a, jusqu’a présent, pu étre trouvé d’interprétation autre que pure- 
ment formelle. 

Le point de départ de cette conception unitaire est dans 1’interpré- 
tation statistique de 1’allométrie de taille convenablement précisée. 


Soit done une population homogéne d’individus ayant atteint leur 
taille définitive et sur chacun desquels peut-étre pratiquée avec pré- 
cision la mesure d’un certain nombre de grandeurs z, y, z.... 
Portons notre attention sur deux de ces variants z et y. Le résultat 
des mesures est un nuage de points (z, y) dont la forme est d’autant 
mieux définie que le nombre d’individus mesurés est plus grand. Nous 
nous limiterons aux cas, heureusement nombreux dans la pratique, ot 
la corrélation entre les deux grandeurs considérées est forte; ce sont 
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d’ailleurs les seuls pour qui la notion d’allométrie puisse présenter un 
véritable intérét. Le nuage figuratif de la population des mesures est 
alors trés allongé. Il s’agit maintenant de préciser sa forme. Sur ce 
point, d’ailleurs capital, il n’est malheureusement pas possible d’arriver 
a une certitude absolue. I] se trouve que, dans bien des eas, le nuage 
peut étre considéré comme sensiblement normal, en ce sens que ses 
points se répartissent 4 peu prés suivant une loi de Laplace-Gauss a 
deux variables et qu’il est allongé suivant une droite. Mais il se trouve 
aussi que, bien souvent, les mémes mesures traduites dans le systéme 
log x, log y donnent encore un nuage de points sensiblement normal, 
allongé suivant une droite. A titre d’exemple de ee fait trés zénéral, 
je citerai les Maia squinado miles adultes, chez qui le coefficient de 
corrélation entre la longueur du propodite de la pince et la longueur 
du céphalothorax est 0.957, et le coefficient de corrélation entre les 
logarithmes de ces grandeurs 0.965, le nuage logarithmique étant peut- 
étre un peu mieux allongé suivant une droite que le nuage arithmétique 
et pouvant étre, au méme titre que ce dernier, considéré comme normal 
L’ambiguité de ce résultat ne tient évidemment qu’a |’insuffisance du 
nombre des individus mesurés et, en augmentant l’importance de 
V’échantillon de la population étudiée, on doit pouvoir faire le choix 
qui s’avére impossible sur un lot de 300 animaux, 4 moins que, ce 
faisant, on ne tombe sur un résultat incompatible avec l’une et |’autre 
des deux hypothéses envisagées. Dans la pratique, avec les populations 
de mesures dont on dispose, |’ineertitude est trés fréquente. 

Il peut certainement arriver que le nuage de points soit plus allongé 
et plus rectiligne en coordonnées arithmétiques qu’en coordonnées 
logarithmiques et que la loi de distribution soit nettement du type 
Laplace-Gauss. II est certain, en revanche, que, dans beaucoup de eas, 
notamment lorsqu’une des mesures porte sur un poids et |’autre sur une 
longueur, le nuage est plus normal lorsqu’il est construit 4 partir du 
logarithme des mesures que lorsqu’il est fait 4 partir de ces mesures 
mémes, la distribution étant alors du type Galton-MacAlister. Mais 
il arrive aussi, et je 1’ai constaté notamment dans 1|’étude des relations 
entre le poids des élytres et celui des mandibules chez Lucanus cervus, 
que le nuage ne soit normal, ni dans un systéme, ni dans l’autre, et ne 
s’allonge suivant une droite que dans sa partie médiane. 

Laissant de cété ces distributions complexes, encore trés insuffisam- 
ment connues, il apparait que la situation est la suivante: dans le plus 
grand nombre des cas, parmi ceux dont |’étude a été la plus compléte 
et la plus rigoureuse, deux lois semblent pouvoir actuellement étre 
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employées aussi légitimement l’une que l’autre. Pour pouvoir justi- 
fier, par des critéres statistiques rigoureux, les choix que 1’on est amené 
a faire entre elles, il faudrait accroitre démesurément |’importance du 
matériel étudié. Encore n’est-il pas certain que ce faisant, on ne 
découvrirait pas qu’une troisiéme loi, plus complexe que les deux 
autres et intermédiaire en quelque sorte entre elles, ne conviendrait pas 
mieux aux faits observés. I] apparait, dans ces conditions, que 1’on est 
actuellement en droit d’exprimer la relation stochastique entre les deux 
variables par une régression logarithmique aussi bien que par une 
régression linéaire, constatation qui légitime en une certaine mesure, 
l’emploi de la relation d’allométrie sans d’ailleurs justifier la place 
qu’a prise cette relation dans la biométrie au détriment de la relation 
linéaire, qui peut cependant prétendre, aussi bien qu’elle, 4 représenter 
les faits observés. Bien que le fait que, dans un certain nombre de cas 
bien établis, la régression logarithmique donne un meilleur ajustement 
que la régression linéaire plaide en faveur de la relation d’allométrie, 
il reste 4 justifier par des arguments biologiques cette préférence 
générale. Nous les trouverons en approfondissant la notion de varia- 
bilité, telle qu’elle est donnée par |’observation de la nature d’une part, 
par la description statistique des faits d’autre part. 

Nous rappelerons tout d’abord qu’aux yeux du biologiste la varia- 
tion relative est plus significative que la variation absolue. Deux or- 
ganes de taille différente d’un méme animal sont considérés comme 
également variables, non pas si leurs dimensions présentent des fluctua- 
tions égales, mais si ces fluctuations ont la méme importance relative ; 
la variabilité s’évalue traditionnellement non pas en unités métriques 
mais en pourcentages. Il apparait naturel, dans ces conditions, de 
dire que deux organes sont également variables si leurs coefficients de 
variabilité sont égaux et de mesurer par ce paramétre l’importance de 
cette variabilité. Il s’agit la d’une définition dont il est naturellement 
impossible de démontrer le bien fondé mais qui apparait en parfait 
accord avec 1’idée que les biologistes se sont fait de la variabilité. Elle 
admet que, dans des organes également variables et de dimensions 
moyennes différentes, l’écart type est proportionnel a la taille, fait en 
parfait accord avec cette constatation empirique que les coefficients de 
variabilité d’organes homologues d’animaux de taille trés différente ne 
présentent pas de variations systématiques avec cette derniére et que 
les objets que le biologiste est porté 4 considérer comme également 
variables ont toujours sensiblement méme coefficient de variabilité. 

Cette définition de la variabilité a plusieurs conséquences intér- 


4 
. 
4 
i 
% 
4 


LA RELATION D’ALLOMETRIE 


23 


essantes. La premiére apparait dés que l’on examine les raisons qui 
peuvent étre données poue expliquer que les variations quantitatives 
d’un étre vivant puissent étre traduites par une loi normale ou sensible- 
ment normale. Ces raisons sont calquées sur celles qui sont utilisées 
pour expliquer la validité de la méme loi dans d’autres domaines et plus 
particuliérement dans ia théorie des erreurs, et il n’y a pas lieu de nous 
y attarder bien longuement. Si des animaux adultes de méme espéce 
et de méme sexe ne sont pas identiques, c’est, d’une part, que leur con- 
stitution génétique présente des différences insoupeonnées et, d’autre 
part, qu’ils n’ont jamais véeu dans des conditions rigoureusement 
identiques. Il peut arriver que les différences de taille soient dues, 
pour la plus grande part, 4 un ou a un petit nombre de facteurs pré- 
dominants, certains génes par exemple; mais on doit admettre qu’en 
général elles résultent de 1’interférence d’innombrables actions, ayant 
chacune un effet trés minime, qui se sont produites 4 un stade on a un 
autre du développement, sous ]’influence des conditions du moment. 
Si nous admettons que, pour la plupart au moins, ces causes de varia- 
tion sont indépendantes et que leurs effets sont additifs, hypothése que 
nous aurons 4 e .miner de plus prés dans un instant, il est possible, 
comme on sait, de démontrer que la distribution finale sera d’autant 
plus proche d’une distribution normale que le nombre des causes de 
variation sera plus grand et leur effet individuel plus petit. La 
démonstration s’étend au eas de distributions 4 deux variables, la 
corrélation s’expliquant, fait pour nous trés important, par le fait 
qu’un grand nombre des causes qui font varier l’une des grandeurs 
agissent également sur |’autre et dans de méme sens. Pour expliquer 
que certaines distributions s’éloignent nettement de la normale il 
suffit d’admettre que certains facteurs préponderants du développe- 
ment n’aient pas une distribution normale dans la population, que, par 
exemple, une partie des animaux ait été soumise, 4 un certain moment, 
a des conditions trés differentes de celles ot restaient placés les autres, 
ou encore qu’une partie des individus soit porteuse d’un géne modifiant 
les dimensions du corps tout entier ou celles de l’une de ses parties.” 
Considérons alors deux races de taille moyenne inégale et supposons 
les si voisines l’une de l’autre dans leur forme et leur structure, que 
jes seules différences existant apparemment entre elles soient d’ordre 
métrique; supposons les en outre également variables, au sens précis 
du terme, l’écart type étant proportionnel a la taille moyenne. D/’aprés 


?Dans ce cas on pourrait s’attendre 4 trouver des distributions de Pearson ou de 
Gram-Charlier. 
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ce qui précéde, il faudra pour qu’il en soit ainsi que le nombre des 
causes perturbatrices, ou la grandeur de leur effet moyen, ou encore le 
nombre et la grandeur des perturbations élémentaires, augmentent avec 
la taille. Plus précisément, il faudra, dans les deux cas extrémes, ou 
que, les perturbations conservant en moyenne la méme importance, leur 
nombre soit sensiblement proportionnel a la racine carrée de la taille, 
ou que, ce nombre restant du méme ordre, |’importance moyenne de 
chacune d’elles doit sensiblement proportionnelle 4 la taille. Il est 
clair que cette deuxiéme hypothése est beaucoup plus raisonnable que 
la premiére et que l’on soit admettre que, lorsque la distribution est 
normale, les causes qui tendent 4 écarter un individu de la moyenne ont 
des effets sensiblement proportionnels 4 cette moyenne. Cette hy- 
pothése dictée par des considérations théoriques a d’ailleurs un fonde- 
ment biologique solide; elle signifie que les facteurs modificateurs de 
la croissance ont des effets multiplicatifs et non pas additifs, qu’ils 
agissent non pas en accroissant ou en diminuant la taille d’une certaine 
fraction de centimétre ou de millimétre, mais en 1’accroissant ou en la 
diminuant d’une certaine fraction de cette taille. Ce fait parait de 
plus en plus solidement établi pour les facteurs génétiques de la crois- 
sance et, en ce qui concerne les facteurs externes, toutes les interpréta- 
tions physiologiques des lois de la croissance globale le postulent plus 
ou moins expressement. 

Autant de raisons pour accepter cette hypothése qu’il devient fort 
intéressant de pousser plus loin. Si |’effet d’un facteur apparait 
comme multiplicatif lorsqu’on compare deux races de taille inégale 
d’une méme espéce, on ne voit pas pourquoi il ne le serait pas aussi a 
l’intérieur d’une méme race lorsque |’on compare les individus les 
plus grands et les plus petits, on ne comprend pas pourquoi |’effet d’un 
facteur sur un individu déterminé serait proportionnel 4 la taille 
moyenne de |’espéce et non pas 4 sa taille propre. On est ainsi conduit, 
par une généralisation toute naturelle, 4 interpréter la variabilité d’une 
population adulte par |’interférence d’innombrables actions élémen- 
taires dont chacune augmente ou diminue d’un pourcentage trés min- 
ime la taille de l’individu sur lequel elle s’exerce. Admettre qu’il en 
est ainsi revient a dire que la distribution des logarithmes des tailles 
doit étre normale, ou que la distribution des tailles elles mémes est du 
type de Galton-MacAlister. Les raisonnements précédents s’étendant 
immédiatement aux distributions 4 plusieurs variables, la régression 
logarithmique apparait comme mieux fondée biologiquement que la 
régression linéaire. 
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Cet argument en faveur de la relation d’allométrie peut étre 
présenté d’une autre facon qui a l’avantage de comporter moins d’hy- 
pothéses. On sait qu’une caractéristique des distributions normales 
a deux variables est d’étre homoscédastiques, e’est-a-dire d’avoir, pour 
chaque variable, un écart type lié indépendant de la valeur particuliére 
de l’autre. Autrement dit, si les individus sont classés selon les valeurs 
croissantes de z, l’écart type de y sera le méme pour les individus les 
plus grands que pour les petits et, de méme, si le classement est fait par 
rapport a y, l’écart type de z restera constant tout au long de la dis- 
tribution. Si la distribution normale est celle des mesures directes, 
il s’ensuit que le coefficient de variabilité des grands individus sera 
inférieur 4 celui des petits. Si la distribution normale est celle des 
logarithmes des longueurs le coefficient de variabilité sera au contraire 
pratiquement indépendant de la taille. Cette remarque ne fournit pas 
pas moins un réel intérét. On n’a jamais remarqué en effet que, dans 
chaque cas particulier, les individus de taille extréme étant gén¢rale- 
ment trop peu nombreux pour qu’un coefficient de variabilité puisse 
étre déterminé sur eux avee quelque sécurité, mais elle n’en présente 
pas moins un réel intérét. On n’a jamais remarqué en effet que, dans 
chaque espéce, les animaux !es plus grands soient moins variables que 
les plus petits; tout porte a croire, au contraire, qu’il n’en est rien et 
cet argument en faveur de la relation d’allométrie mérite d’étre 
retenu.* 

Tenu compte de tout ce qui précéde, nous opterons pour la relation 
d’allométrie contre la relation linéaire, en sachant bien cependant que 
ce choix raisonnable n’est pas obligatoire, puisqu’en matiére de liaison 
stochastique, il reste toujours place pour un large arbitraire dans le 
choix des formes analytiques. Nous pourrons méme dans certains cas 
user de cette latitude pour utiliser alternativement, selon le but pour- 
suivi, l’une et l’autre des deux relations coneurrentes.* 


*On devrait constater par exemple que le coefficient de variabilité du propodite de la 
pince des Maia squinado males est de 9.8 pour les animaux mesurant 140 mm et de 4.1 
pour ceux qui atteignent 220mm, avec une valeur voisine de 5.8 pour les animaux de 
taille moyenne 180mm, les changements de valeur de V dans un intervalle d’environ 2 
écarts types de part et d’autre de la moyenne etant particuli¢rment importants dans ce cas 
ou la variabilité de lorgane considéré est trés grande (V = 21). Pour le méropodite du 
4éme péréiopode dont la variabilité, encore, forte, est cependant plus normale (V = 12.7) 
les coefficients de variabilité correspondant aux mémes dimensions du céphalothorax 
devraient étre respectivement 4.7 ; 2.8 et 3.5. Rien de tel n’a été observé. 

‘Voir en particulier l'usage que j’ai fait de cette possibilité dans une comparaison 
biométrique de deux espéces du genre Maia, M. squinado et M. verrucosa (Comptes 
Rendus de VAcademie des Sciences, T 204, p. 67, 1937). 


| 
= 
| 
| 
> 
‘| 


26 BIOMETRICS - MARCH 1948 


Le seul argument qui pourrait incliner 4 préférer la traduction de 
la relation existant entre les deux variants par une formule linéaire 
plutét que par une formule logarithmique est que, si deux parties d’un 
méme organe comparées 4 une méme grandeur de référence obéissent 
respectivement a des relations d’allométrie de constantes a; b,; a2 be, il 
est impossible, si a;, et a, sont différents, que l’organe tout entier obéisse 
a une relation de la méme forme. La difficulté ne se présente pas pour 
la relation linéaire ov il suffit de poser : 


Y=9i+ @=0,+4, 


L’objection serait valable s’il s’agissait de relations fonctionnelles 
puisque dans ce cas, les seules relations répondant 4 la condition 
d’additivité doivent étre de la forme y = a¢ (x) +b, condition que ne 
remplit pas la relation d’allométrie, ni aucune de celles que 1’on pour- 
rait songer 4 lui opposer, puisque la réciprocité entre y et xz exigerait 
que l’on ait en méme temps z= a’ ¢ (y) +b’, ce qui nous fait retomber 
sur la relation linéaire. J’ai montré, il y a plus de quinze ans, que 
cette objection est sans valeur dans le probléme qui nous occupe, ot la 
relation cherchée est stochastique. Sauf en des cas exceptionnels, on 
obtiendra une relation satisfaisante en prenant pour a la moyenne 
pondérée de a, et a2: 


ort (a; — az) (91 — 
Yi+ Yeo 2 Yit Yeo 


La différence entre cette valeur moyenne de a et celle qui conviendrait 
pour l’extrémité de la distribution ot 1’écart est le plus grand est 
voisine de V, (a,—a,)*. Elle est presque toujours dans la pratique 
inférieure 4 l’erreur que l’on peut commettre dans la détermination 
d’une constante d’allométrie, qui est, nous le verrons, de l’ordre de 


a\/1-r/\JN. 


II 


Avant de poursuivre notre analyse de la notion d’allométrie, il nous 
faut traiter un probléme fort important en lui-méme et dont la solu- 
tion est moins simple qu’on ne l’a cru. Une fois admis qu’une rela- 
tion d’allométrie peut exister entre deux variables, il faut, pour véri- 
fier cette hypothése, disposer d’un méthode réguliére qui permette de 
déterminer les constantes qui la définissent et la précision des chiffres 
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obtenus.’ Une premiére question se pose immédiatement. L’ajuste- 
ment doit-il étre pratiqué sur l’expression y = bz* ou sur |’expression 
équivalente log y = a log z+log b. Les arguments donnés en faveur de 
l’ajustement de la premiére formule ne sont pas particuliérement con- 
vaineants. Les ecaleuls sont d’une conduite plus difficile et donnent 
d’ordinaire, au surplus, des résultats fort peu différents de ceux que 
l’on obtient en utilisant la relation sous sa forme logarithmique. 
Comme par ailleurs nous avons conclu de 1’étude précédente que la 
distribution de Y = log y et X = loz x doit en principe étre normale, il 
apparait naturel de faire tous les calculs sur la deuxiéme expression 
que nous écrirons Y = a X + B et dans laquelle il s’agit de déterminer a 
et B. 
La premiére solution qui se présente a l’esprit est celle qui consiste 
a identifier la droite cherchée avec la ligne de régression de Y en X. 
Dans le cas of nous nous sommes plaecés par hypothése, celui de 1’étude 
d’un échantillon représentatif de la population globale, cette solution 
ne peut pourtant pas étre retenue. Elle suppose en effet qu’on est en 
droit de faire jouer un réle différent aux deux variables, l’une X étant 
considérée comme indépendante et l’autre Y comme dépendante. Il 
peut évidemment arriver que cette hypothése soit acceptable, si, par 
exemple Y est la mesure d’un organe de petite taille et X celle du corps 
tout entier, mais en général cependant, X et Y jouent des réles symétri- 
ques. Sil’on doit comparer deux mesures linéaires du corps d’un ani- 
mal, longueur et largeur du céphalothorax d’un crabe par exemple, ou 
de deux appendices homologues, on ne voit pas bien de quel droit on 
considérerait l’une plutét que l’autre comme représentant la variable 
indépendante, ni par conséquent comment on pourrait choisir entre les 
deux lignes de régression. L’emploi de ces derniéres ne se justifie que 
si, pour des raisons pratiques, on veut pouvoir calculer 4 partir de l’une 
des dimensions prise, par convention, comme référence, la valeur la 
plus probable de l’autre, mais dans ce cas on ne travaillera sans doute 
pas sur un échantillon représentatif de la population et 1’on ne se trou- 
vera done pas dans le cas considéré actuellement. Cet emploi n’est 
généralement pas licite, les deux variables étant égales en dignité, a 
moins que l’on ne se résigne 4 donner au probléme posé, non pas une, 


5 Le probléme a été étudié en détail par Kavanagh et Richards (Proceedings of the 
Rochester Academy of Science, 150, 1942) qui ont montré que le traitement mathéma- 
tique doit varier avec la nature des données de l’observation et le but poursuivi. Je ren- 
verrai 4 cet excellent travail pour tout ce qui touche A Vhistorique de la question, a 
Vinterprétation de 6 et & sa prétendue corrélation avec a. Ils ont justement insisté sur 
la nécessité qu’il y a dans certains cas, A faire jouer A y et z un réle symétrique dans les 
calculs, mais sont restés dans le cadre des techniques traditionnelles. 
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mais deux solutions qu’il faudra toujours considérer simultanément. 
Si nous voulons arriver 4 une droite unique, la ‘‘meilleure’’ possible, 
nous devons renoncer au critére usuel de la qualité d’un ajustement et 
a l’emploi de la méthode des moindres carrés dans sa forme tradition- 
nelle, rien n’autorisant 4 rendre minimum plutét la somme des carrés 
des écarts des Y 4 X constant que celle des X & Y constant. Le pro- 
bléme est trés voisin de celui qui se pose lorsqu’il s’agit d’ajuster une 
droite 4 un ensemble de couples de mesures ou les deux variables sont 
sujettes 4 l’erreur, probléme qui n’a été que peu travaillé jusqu’a 
une date récente et qui n’a pas encore regu de solution vraiment satis- 
faisante. La meilleure parait encore celle qui rend minimum la somme 
des carrés des distances des points a la droite; la droite obtenue se con- 
fond, on le sait, avec le grand axe des ellipses d’égale distribution. Du 
point de vue géométrique la solution est parfaitement raisonnable; du 
point de vue statistique elle 1’est moins, car on ne voit pas trés bien ce 
que peut représenter en terme d’erreur la distance d’un point figuratif 
aladroite. Il m’asemblé qu’un autre principe de minimum pourrait 
étre utilement envisagé et je propose de représenter la relation d’allo- 
métrie par la droite (D) qui rend minimum la somme du produit des 
écarts de Y pour X constant et de X pour Y constant, ce qui revient a 
rendre minimum la somme des aires des triangles rectangles ayant pour 
hypothénuse commune la droite cherchée, deux cédtés paralléles aux 
axes, et leurs sommets aux points figuratifs. 

Un tel procédé n’est au fond que la transposition dans le domaine 
du caleul d’une technique d’interpolation graphique employée sys- 
tématiquement par certains dessinateurs et que l’on utilise inconsciem- 
ment lorsque l’on essaie de placer au mieux une droite au travers d’un 
nuage de points approximativement alignés. 

On voit sans peine que la droite cherchée (D) passe par le centre de 
gravité X, Y, de la distribution. Sa pente est donnée par la valeur 
de a qui rend minimum la somme des produits [(Y- Y) -a(X-X)] 
[(X -X) -(¥Y-Y)/a], e’est a dire par a? = cy?/ox*. La solution posi- 
tive, a = cy/ox donne la paramétre cherché, la solution négative, a’ = —a, 


*Un autre inconvénient de la solution classique est que la droite obtenue dépend des 
unités adoptées pour la mesure des deux grandeurs comparées I] suffit d’ailleurs pour 
y parer de normaliser les coordonnées, c’est A dire de prendre pour unités de mesure les 
écarts types des deux grandeurs. Mais si un tel artifice peut, A la rigueur, se défendre 
lorsque la comparaison porte sur les mesures mémes, il ne saurait raisonnablement étre 
utilisé lorsque les variables intervenant dans les calculs sont des logarithmes. Nous 
verrons cependant que la droit de moindre distance ainsi obtenue se confond avec celle & 
laquelle nous arrivons par des considérations trés différentes et beaucoup plus acceptables. 
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Fig. 3 
Définitions comparées de la droite d’allométrie D et des deux droites de re- 
gression B et C et des points correspondant & X, Y, sur D (X’, Y’,) et sur les 
droites de régression (X, Y”, et X”, Y,). 


est la pente de la droite (D’) qui rend maximum cette somme de pro- 
duits. 
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L’équation de la droite (D) s’écrira: 


— 
ox ox 


mais pourra aussi se mettre sous la forme symétrique: 


Y-Y_ X-X 
oy ox 


La variance de la constante d’équilibre a est ca* = a?(1—r*)/N, r étant 
le coefficient de corrélation de X et de Y; elle est égale a celle du co- 
efficient de régression de Y en X. Les droites (D) et (D’) constituent 
un couple de directions conjuguées pour les ellipses d’égale distribu- 
tion et le point X’,, Y’, de la droite (D) correspondant aux mesures 
X;, Y;, relatives 4 un individu particulier est situé sur la paralléle a 
(D’) menée par X,, Y,. Ses coordonnées sont: 


or 


et l’on peut aussi écrire: 


ox oY 2 


La variance des écarts entre les points figuratifs et les points de la 
droite (D) qui leur correspondent est constante et égale & (cx? + cy’) 
(1-1r)/2. La mesure de la surface moyenne des triangles d’ajustement 
joue par ailleurs dans nos calculs le méme réle que la variance liée de 
la méthode classique de la régression et peut étre qualifiée de covariance 
liée. Elle est proportionnelle 4 la grandeur précédente et constam- 
ment égale en tous points de la distribution 4 cy cy (1-r). On voit 
que les formules qui précédent permettent de faire correspondre sans 
ambiguité 4 tout individu réel de la population l’individu normal 
fictif qui s’en rapproche le plus et de mesurer en écarts réduits leur 
degré de dissemblance qui peut étre défini par 1’expression : 
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Fig. 4 
Les droites D et D’ mises en place dans une des ellipses d’égale distribution ; 
B et C, les deux droites de regression et A, 1’axe principal (a = 1,285, r=0,71). 


L’usage qui vient d’étre fait de la formule d’allométrie pour définir 
Vindividu théorique correspondant 4 chacun des individus rels étu- 
diés, est celui pour lequel elle a été établie et auquel elle est le plus 
exactement adaptée. Mais, dans la pratique, elle doit étre utilisée beau- 
coup plus fréquemment 4 estimer la valeur de Y qui correspond nor- 
malement 4 une certaine valeur de X ou inversement a estimer la 
valeur de X correspondant normalement 4 une certaine valeur de Y. 
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Elle remplace ainsi, suivant le cas, 1’une ou |’autre des deux lignes de 
régression de Y en X ou de X en Y. y 
L’estimation de Y qu’elle apporte différe de celle que fournit la 


ligne de régression de Y en X de (1-1) oy X-<X Par ailleurs la 
ox 

variance de ]’écart entre les Y observés et les Y calculés par cette for- 

mule pour une valeur donnée de X est : 


a [143-8 


l+r 
que l’on peut aussi écrire : 


Bert (1-r) (1-r)?| 


Elle dépend donc de (X-X), sa valeur moyenne est 2cy* (1-r). 
Elle n’est égale 4 la variance liée de Y, variance de 1’écart entre les Y 
observés et les Y calculés par la ligne de régression, qui est égale, quel 
que soit X, a cy? (1—r*) qu’a sa valeur minimum atteinte pour X = X. 
L’amplitude de ces variations est d’ailleurs faible. Il est exceptionnel 
que dans les échantillons de dimensions usuelles elle dépasse sensible- 
ment, pour les individus les plus grands ou les plus petits, les erreurs 
commises normalement dans la détermination de la variance liée. Pour 
la méme raison, et bien que les estimations fournies par la droite d’allo- 
métrie présentent des écarts systématiques avec celles que permet la 
droite de régression, il est rare, lorsque la corrélation est forte, comme 
nous l’avons constamment supposé,- que ces écarts soient notablement 
plus amples que ceux qui se rencontrent normalement entre points ob- 
servés et points calculés par l’une ou |’autre droite. Il reste évidem- 
ment possible que, dans certains cas, et pour certains usages, il soit pré- 
férable d’utiliser la ligne de régression classique, mais cette reserve ne 
diminue en rien 1’intérét de la droite que nous avons définie. 

Une des caractéristiques les plus remarquables de cette droite est 
que, parmi toutes celles qui pourraient étre employées 4 décrire la 
liaison de X et de Y, elle est la seule dont la pente soit indépendante de 
la valeur du coefficient de corrélation de X et Y. Cette caractéristique, 
qui n’aurait peut-étre pas grand intérét dans 1’étude d’une distribution 
quelconque, est fort avantageuse dans une étude de morphologie quan- 
titative ot, comme nots le verrons, la valeur numérique du coefficient 
de corrélation ne présente pas d’intérét biologique. I] importe en 
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Fig. 5 
Allométrie de taille chez les Maia Squinados miles adultes. Longueur du 
propodite de la pince en fonction de la longueur du méropodiie du dernier pereiopode 
locomoteur chez 301 individus. Coordo logarithmiques (a = 1,57, r=0,963). 


effet de bien comprendre que la loi que 1’on essaie de dégager de 1’en- 
semble des observations est, en définitive, celle que l’on pourrait raison- 
nablement espérer voir se manifester si la corrélation tendait vers la 
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perfection. Imaginons en effet une suite de distributions & deux varia- 
bles & distributions marginales invariables et & corrélation croissante. 
On verra alors les lignes de régression se rapprocher 1’une de 1’autre et 
tendre, en méme temps que |’axe principal vers une position limite qui 
est précisemment celle de la droite (D) qui, elle, reste fixe.’ C’est 
cette permanence qui, avec sa construction symétrique par rapport aux 
deux variables, lui donne pour nous un intérét particulier, qu’accroit 
encore la facilité de son calcul. 

Il importe cependant de préciser que ce calcul n’est valable qu’au- 
tant que le lot étudié constitue un échantillon représentatif de la popu- 
lation. Il cesserait de pouvoir étre utilisé si les individus mesurés 
avaient été sélectionnés d’aprés les valeurs de l’une des variables. La 
seule méthode recommandable dans ce cas est la méthode des moindres 
carrés classique. 

Un avantage important de la méthode réguliére lorsqu’elle peut étre 
employée et que l’on étudie simultanément plusieurs organes, est de per- 
mettre de donner en une fois toutes les relations d’allométrie qui peu- 
vent exister entre toutes les variables prises deux 4 deux. I] suffit pour 
cela d’écrire les équations sous leur forme symétrique: 


x-X¥ Y-Y¥ 
ox oy on 
Toutes les relations cherchées se trouvent ainsi représentées, si 1’on 
compare p organes,-par une seule droite de l’espace 4 p dimensions, 
droite dont les diverses projections représentent les relations des or- 
ganes pris deux a deux, trois a trois, ete. . . . - 


' Le pente de la droite (D) est intermédiaire entre celle de l’axe principal et celle de 
la bissectrice des lignes de régression. Dés que la corrélation est forte ces trois droites 
sont trés voisines. Avec r=0.90 et a= 0.900 la pente dé l’axe principal est 0.889 et celle 
de la bissectrice 0.9005 ; avec r = 0.95 et a = 0.900, les pentes correspondantes sont 0.895 et 
0.9001. La moyenne des pentes des deux droites de régression donne elle-méme une trés 
bonne approximation respectivement 0.905 et 0.901. Dans un cas comme celui des Maia 
ou toutes les corrélations dépassent 0.95, toutes ces droftes se confondent pratiquement, 
les erreurs d’énchantillonnage pouvant dépasser plusieurs milliémes. En intervertissant le 
réle des variables dans !’exemple précédent, on verrait que les conclusions valables pour 
Yallométrie a= 0.900 le sont également pour l’allométrie inverse a= 1.111. 

I est facile de voir que lorsqu’ll y a isométrie (a =1), axe principal, bissectrices des 
lignes de regression et droite d’allométrie se confondent. 

Nous obtiendrions la méme superposition pour une valeur de a quelconque & la condi- 
tion de faire usage de coordonnées normalisées. On peut done dire, en un certain sens, 
que la solution proposée ici se confond avec celle qu’a préconisée Pearson, sans que la 
justification qu’il en avait donnée ait jamais paru pleinement satisfaisante. Il importe 
de noter & ce propos que le principe de minimum dont nous avons fait usage et qui nous 
parait constituer la meilleure justification de la droite d’allométrie, ne semble avoir été 
invoqué, ni par Pearson, ni par aucun de cenx qui ont étudié plus recemment le méme 
probleme. 
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Sans vouloir développer ici cette extension des résultats précédem- 
ment obtenus, il est intéressant de noter que, si l’on prend en considéra- 
tion l’ensemble des p organes étudiés, on peut, comme nous |’avons 
fait pour deux organes seulement, faire correspondre 4 chaque indi- 
vidu réel, X, Y, Z,.. . un individu normal caractérisé par des mesures 
X’,, Y's, . . . telles que: 


Y.-Y¥_2,-Z Z,-Z | 


ox oy oz p 


Tous les calculs précédents ont été faits sur les logarithmes et cette 
procédure est le seule qui soit reeommandable. Des formules connues 
permettent d’en déduire les constantes qui caractérisent la distribution 
des mesures elles-mémes. Dans la pratique les formules ci-dessous 
donnent une approximation toujours suffisante. 

En appelant #’, le nombre dont le logarithme est égal 4 la valeur 
armoyenne de log y, on a approximativement : 


=e" (1+cy? /2) 
oy=Yoy (1+ 30,7/4) V, = oy (1+ cy?/4) 


oy ne peut guére dépasser 0.2, ce qui correspond déja 4 un rapport de 1 
4 4 entre les dimensions du variant pour les individus de taille extréme. 
On voit que l’on peut écrire 4 moins d’un centiéme prés : 


Vy=oy 


approximation toujours suffisante puisque |’erreur d’énchantillonnage 
commise sur le détermination de 1’un ou 1’autre de ces paramétres peut 
dépasser largement le centiéme tant que 1’effectif de 1’échantillon est 
inférieur 4 quelques milliers. On aurait inversement, si les caleuls 
avaient été faits sur les mesures mémes : 


Y=logg-V,7/2 oy = V,(1-V,?/2) 
ou pratiquement : 


oy=Vy 


Comme nous devions nous y attendre, la constante d’équilibre est le 
quotient des coefficients de variabilité des deux grandeurs comparées. 
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Nous pouvons maintenant reprendre notre étude du probléme bio- 
logique de 1’allométrie. 

Nous avons vu, dés le début de notre exposé, que l’allométrie de 
taille était trés étroitement liée 4 la variabilité et nous venons de ra- 
mener & un seul ces deux problémes. Définir la constante d’equilibre 
« par le quotient des deux coefficients de variabilité revient 4 dire 
qu’il y a allométrie dés que les grandeurs que 1’on compare sont in- 
également variables et que l’isométrie ne ‘peut exister que si les deux 
variabilités sont les mémes. Nous avons interprété cette variabilité 
d’une population adulte par le jeu d’innombrables causes dont chacune 
@ pour effet d’augmenter ou de diminuer d’un pourcentage trés 
minime les dimensions d’un organe ou d’un autre. On sait, par ail- 
leurs, et nous allons avoir |’occasion de le préciser, qu’il y aura corré- 
lation entre deux organes si certains facteurs agissent simultanément 
sur l’un et sur l’autre, la corrélation étant d’autant plus compléte que 
le nombre des facteurs est plus grand et leur action individuelle plus 
importante. Ii est clair enfin que |’allométrie ne serait parfaite que 
si la corrélation |’était aussi, les variances liées ne devenant nulles qu’a 
cette occasion. Plus généralement, la relation d’allométrie sera d’au- 
tant plus exactement suivie que la corrélation entre les deux grandeurs 
comparées sera plus rigide et qu’il y aura moins de causes faisant varier 
un des organes indépendamment de |’autre. 

Il serait par conséquent tout a fait déraisonnable de chercher une 
explication de |’allométrie de deux organes dans le jeu de certains fac- 
teurs qui agiraient sur l’un des organes et non pas sur |’autre, |’inter- 
vention de tels mécanismes ne pouvant, au contraire, que troubler la 
pureté du phénoméne. S’il existe des facteurs de croissance agissant 
plus ou moins électivement sur telle ou telle partie de corps, leur action 
ne devra étre invoquée que pour justifier ]’existence d’écarts plus ou 
moins importants entre les prévisions permises par la relation d’allo- 
métrie et la réalité, l’allométrie méme s’expliquant entiérement par 
l’inégale sensibilité des deux variants au méme ensemble de facteurs. 
Si, par exemple, la longueur de propodite de la pince des Maia males 
comparée a la longueur du céphalothorax a une constante d’équilibre 
de 1.93, c’est qu’en moyenne le propodite de la pince des Maia est 1.93 
fois plus sensible que le céphalothorax & toutes les actions qui peuvent 
s’exercer au cours du développement en le favorisant, ou le défavori- 
sant. 

Il est facile de préciser quantitativement ces indications en raison- 
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nant sur un modéle extrémement simplifié, mais cependant fort in- 
structif. Imaginons le cas idéalement simple of deux organes com- 
parés réagiraient aux mémes causes indépendantes de variation et & 
aucune autre, et ol, pour chacune d’elles, 1’effet serait proportionnel & 
l’intensité, le facteur de proportionnalité étant le méme pour tous les 
facteurs et différent pour les deux organes. Adoptons pour chacun 
des facteurs un systéme de mesures tel que sa valeur moyenne soit nulle. 
Nous pourrons alors écrire, U2, . . . , étant les facteurs de varia- 
tion supposés trés nombreux : 


X =X+ du, + Ay t+. Ally 


Nous savons que ces deux distributions sont normales ou approxima- 
tivement normales et nous avons: 


= (os,? + os,” + os,” 


La constante d’équilibre sera alors a= »/A, quotient des facteurs de 
proportionnalité. Le coefficient de corrélation sera égal a l’unité cur: 


= 


Og Vy 


Dans ce cas schématique, allométrie et corrélation seront simul- 
tanément parfaites, parce que toute la variabilité est conditionnée dans 
les deux organes par les mémes causes. 

Imaginons maintenant qu’on changement dans les conditions 
d’éxistence de l’espéce, ou un changement génétique, fasse disparaitre 
un groupe de causes de variabilité, mais leur en substitue deux autres, 
indépendants, qui exerceraient des actions de méme nature et de méme 
amplitude que les facteurs disparus, le premier agissant électivement 
sur l’un des organes, le second agissant électivement sur l’autre. Par 
hypothése, les variances de X et Y restent les mémes, mais le coefficient 
de corrélation diminuera, les nouveaux facteurs n’intervenant pas dans 
son numérateur, tandis que son dénominateur reste inchangé. Si 
U1, Ua, Us... , Wi, We, Ws . . . sont ces facteurs, qui exercent re- 
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spectivement, par hypothése, des actions Av,, Av2, Avs . . . SUT ¥ et 
pW, pWe, pw, ... sur Y, on trouve: 


et: 
ox’ 


On vérifie aisément que les numérateurs de ces deux derniéres 
expressions représentent les valeurs moyennes des variances liées de 
X et de Y ce qui nous redonne les expressions connues de ces grandeurs, 
2(1-1r)ox? et 2(1-1) cy’. 

On voit que, dans cet exemple, allométrie et corrélation sont indé- 
pendantes l’une de l’autre. L’allométrie reste constante parce que, par 
hypothése, l’effet résultant de l’ensemble des facteurs reste le méme, 
lorsque certains d’entre eux sont remplacés par d’autres; le coefficient 
de corrélation change avec le nombre des facteurs communs aux deux 
organes. Lorsque le nombre des facteurs communs diminue, la dis- 
persion des points figuratifs autour de la droite d’allométrie augmente, 
ce que traduit l’accroissement des variances liées de X et de Y; la 
détermination de a devient en méme temps plus imprécise, comme le 
montre l]’expression de sa variance, a? (1-r*)/N. La valeur de la 
constante d’equilibre a est ici égale & w/A, c’est a dire au quotient des 
deux coefficients qui mesurent la sensibilité respective des deux organes 
aux actions perturbatrices. Nous voyons ainsi apparaitre, avec une 
netteté toute particuliére sur ce schéma simplifié a 1’extréme, la carac- 
téristique essentielle du phénoméne d’allométrie que nous indiquions 
tout l’heure propos d’un cas concret : l’allométrie de taille exprime 
l’inégale variabilite de deux organes, témoignage elle-méme de leur 
inégale sensibilité 4 1’ensemble des facteurs perturbateurs de la crois- 
sance. 

Nous pourrions aisément compliquer le schéma précédent et nous 
rapprocher ainsi progressivement de la réalité. On congoit sans qu’il 
soit besoin d’explication, qu’il peut arriver, qu’il doit méme arriver, 
que les différent facteurs aient une importance inégale, que chacun 
d’eux ait ses coefficients d’action propres. Les résultats changeront 
naturellement, mais l’essentiel en subsistera. L’allométrie restera 
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imputable a 1’ensemble des facteurs de la variabilité, la corrélation ne 
dépendra toujours que des seuls facteurs communs aux deux organes 
comparés ; la variabilité résiduelle, que mesure la variance liée et qui 
est responsable des écarts des points observés et de la droite figura- 
tive, restera imputable aux facteurs qui agissent sur 1’un ou sur 1’autre 
des deux organes, mais non pas sur les deux. 

Si notre étude se limite 4 la comparaison de deux organes, nous ne 
pouvons aller plus loin; mais si nous étudions plusieurs organes du 
méme animal nous pouvons constater que chaque organe a sa varia- 
bilité propre et que, suivant les couples d’organes comparés, la correla- 
tion est plus ou moins forte. On peut alors se demander si, en dehors 
des facteurs généraux qui agiraient sur toutes les parties de l’organisme 
et expliqueraient qu’entre deux organes quelconques existe toujours 
un certain degré de corrélation, n’existeraient pas d’autres facteurs 
agissant de facon plus ou moins élective sur la croissance de telle ou 
telle partie de l’animal. Ceux-ci pourraient d’ailleurs appartenir a 
plusieurs catégories, les uns étendant leur action 4 tout un systéme 
d’organes, les autres n’exercant leur influence que sur une partie plus 
limitée du corps. Une telle étude, sans apporter aucun élément 
nouveau 4 la connaissance de |’allométrie méme, permettrait de décéler 
certaines influences qui l’empéchent de se manifester dans tous les cas 
avec la méme rigueur. On peut espérer d’ailleurs qu’elle pourrait 
apporter quelque lumiére sur des faits qui n’auraient méme pas été 
soupconnés sans elle.*® 

Posé dans ces termes, le probléme est typiquement du ressort de 
l’analyse factorielle. Cette technique statistique, actuellement d’un 


emploi constant en psychologie, n’a encore été utilisée que trés rare- - 


ment en biométrie bien que, dés 1932, S. Wright en ait montré ]’in- 
térét. Il faut espérer que les méthodes de calcul mises au point récem- 
ment par Delaporte en généraliseront 1’emploi. J’en ai moi-méme 
fait usage dans un travail consacré 4 1’étude des variants sexuels d’un 
Crustacé et j’indiquerai ci-aprés les principes qui m’ont guidé dans 
cette étude. Ils ne sont d’ailleurs que la généralisation de ceux que 
nous avons déja mis en oeuvre. 


8 Je citerai A titre d’exemple, le probléme de la taille cellulaire. Dans un travail cons- 
sacré & l'étude des facteurs qui la conditionnent chez les Mammiféres (C. R. Soc. Biol., 
T. 135, 1941, p. 662, 750 et 1309), j’ai pu montrer en premier lieu que, d’une facon trés 
générale, les dimensions du noyau sont unies a celles de la cellule par une relation de la 
forme N = bC® ou b est une constante caractéristique du type cellulaire et ou a est voisin 
de 0.6 ou 0.7. J'ai établi ensuite que la taille moyenne d’une cellule dépend de deux fac- 
teurs, l’un général qui a des valeurs approximativement égales pour des animaux de 
méme taille, quelle que soit leur espéce, l’autre spécial, qui ne a que de la nature de 
Vanimal et non pas de sa grandeur. 
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Imaginons que nous ayons étudié un certain nombre d’organes et 
que les facteurs qui agissent sur eux se répartissent en un certain nom- 
bre de groupes indépendants entre eux U2, .. . U's, . . 
w”’,,.... Pour plus de commodité nous désignerons chaque groupe 
par une seule lettre, de sorte que nous pourrons écrire pour chaque 
organe : 


un ou plusieurs des coefficients A,, A’1, Ax, . . . pouvant étre nuls. 

Le nombre des groupes de facteurs n’est pas défini a priori et peut 

varier avec le nombre d’organes étudiés, deux groupes pouvant agir 

simultanément sur plusieurs organes et séparément sur certains autres. 
On aura alors: 


oy,” = ov’ + r’,* Cy? + why 


la fraction de la variance totale attribuable a U est 7? = A,’oy?/cy,* et 11 
représente le coefficient de corrélation r,y de Y, et de U. On montre 
sans difficulté que le coefficient de corrélation entre deux organes Y, et 
Y, est: 

To9 = Tp0T gu + + - 


Les coefficients de corrélation r,, pouvant étre mesurés directement 
il s’agit de calculer 4 partir d’eux les coefficients . . Tu, 
fv Il n’existe pas de méthode générale de résolution d’un tel systéme 
d’équations et l’on doit procéder 4 une série de tatonnements méthod- . 
iques, en essayant successivement des schémas a 1, 2, 3 . . facteurs, 
jusqu’a ce que 1’on arrive a retrouver.les rpy avec les pg. 

L’application de cette méthode® m’a permis de reconnaitre chez les 
Masa un facteur général qui conditionne la taille globale, évaluée par les 
dimensions du céphalothorax, dont la mesure L est la meilleuer esti- 
mation et qui contréle 92 4 94% de la variance totale des appendices ; 
un facteur de groupe G qui agit sur l’ensemble des appendices et 
contréle 3 4 6% de leur variance; un facteur regional R qui agit 
électivement sur les appendices antérieurs dont il contréle 3% environ 
de la variance. La variance résiduelle, qui est de l’ordre de 1% de la 
variance totale, est imputable 4 des facteurs locaux et ne peut étre 
évaluée que par différence. 

Il est possible & partir de 1’ensemble des mesures des appendices de 
cealculer les valeurs numériques de G et de F relatives 4 chaque individu, 

* Voir Biotypologie, T. 6, pp. 73-98, 1938. 
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et d’établir, d’autre part, des formules de prévision plus complétes que 
la relation d’allométrie, en ce qu’elles font intervenir simultanement 
L,G,®. On constate alors qu’il est possible de recalculer, 4 partir de 
eet ensemble de trois nombres, les dimensions de chaque individu avec 
une précision supérieure au centiéme, les coefficients de variabilité 
résiduelle étant tous inférieurs au dixiéme des coefficients de variabilité 
initiaux. 

On ne peut pas prévoir, tant que d’autres recherches sur ce probléme 
n’auront pas été faites, quel peutétre l’avenir de ce genre d’étude. 
Il n’est pas interdit d’espérer qu’elles pourront apporter quelques 
lumiéres sur les ow la relation d’allométrie se vérifie imparfaitement, 
soit que le nuage de points figuratifs soit trop étalé, soit qu’il s’incurve 
a l’une de ses extrémités. 

IV 


Il s’agit de montrer maintenant que le probléme de 1’allométrie de 
croissance n’est pas différent au fond de celui de 1’allométrie de taille. 

Nous remarquerons en premier lieu que rien de ce qui a été dit plus 
haut de l’allométrie de taille n’implique que les animaux comparés 
aient définitivement cessé de croitre, mais que nous avons seulement 
supposé que la population étudiée ne comportait que des individus 
ayant atteint le méme stade de développement. Nos conclusions vau- 
draient, par exemple, tout autant pour chacun des stades larvaires d’un 
Insecte que pour son stade imaginal. C’est a partir de cette trés simple 
remarque que nous allons pouvoir construire la démonstration cherchée 
que nous établirons tout d’abord sur le développement d’un Arthro- 
pode, la discontinuité de la croissance dans cet embranchement rendant 
les faits particuliérement clairs. 

Considérons les stades successifs du développement d’un de ces 
animaux, un Crustacé par exemple, en les choisissant dans une méme 
phase de la croissance, c’est 4 dire dans une période de 1’existence ot 
ne se produit aucun remaniement organique, ni aucun changement 
notable dans le fonctionnement physiologique des différents appareils. 
Supposons, plus précisément, que nous sommes dans les conditions ot 
une étude de la croissance relative a un sens, ¢’est & dire que la forme 
d’un individu ne dépend que de sa taille. Il peut arriver que la 
population de chacun de ces stades soit assez peu variable et que le 
taux d’accroissement lors de la mue soit assez grand pour que les 
représentants de deux stades consécutifs aient des dimensions qui ne 
se chevauchent pas. En général cependant les individus les plus grands 
du premier stade ont des dimensions supérieures 4 celles des individus 
les plus petits du second stade. Autrement dit, les nuages figuratifs 
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Fig. 6 
ope Schéma correspondant a celui de la fig. 4, dans l’hypothése ot la population 
Ho comporte en nombre égal les représentants de deux stades consécutifs du dévelop- 
pement. Voir le texte. 
r des deux populations de mesures de x et de y se recouvrent partielle- 
| ment. Comme nous avons supposé expressément que dans l’espéce 
‘i étudiée la forme ne dépend pas directement de 1’Age, les deux séries de 


7 mesures de la zone de recouvrement doivent nécessairement corre- 
fe spondre 4 la méme loi de croissance et donner des resultais aussi con- 
cordants que si elles avaient été pratiquées dans la méme zone sur deux 
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populations appartenant 4 un méme stade. II résulte de 1a que, si une 
courbe puissance traduit dans cette zone pour la premiére population 
la relation existant entre les deux variants, elle doit la traduire tout 
aussi bien pourladeuxiéme. Les relations d’allométrie f correspondant 
a la premiére population et f, correspondant a la deuxiéme, qui coinci- 
dent dans la zone de chevauchement, coincideront alors dans toute leur 
étendue. Il en serait de méme pour tous les stades qui précédent ou 
qui suivent les deux stades considéres jusqu’a présent, a la seule condi- 
tion que dans toute cette étape du developpement, on soit en droit de 
considérer la forme comme fonction de la taille. La loi de croissance 
f sera en ravanche remplacée par une loi différente ¢ lorsque le dévelop- 
pement sera arrivé a un stade tel que, dans la zone de recouvrement, 
les deux populations doivent étre considérées comme différentes. 

Le probléme de la recherche d’une loi de croissance relative valable 
pour une des étapes du développement d’un Arthropode se raméne 
ainsi, dans les hypothéses ot: nous sommes placés, 4 celui de la recherche 
de la relation existant entre deux organes pour les animaux d’un méme 
stade choisi arbitrairement dans cette étape. Ainsi s’explique, de la 
facon la plus naturelle, ]’identité formelle des relations qui permettent 
de comparer les formes successives realisées au cours de développement, 
ou les formes diverses que peuvent présenter suivant leur taille les 
:nimaux d’un méme stade. Ainsi s’explique aussi trés simplement, le 
fait que le passage d’une <tape a la suivante se marque sur la courbe de 
croissanee relative par un point anguleux ou une discontinuité, 
puisqu’il n’y a aucune raison pour que les courbes f et ¢ se raccordent 
tangentiellement l’une a 1’autre. 

L’élément nouveau que nous venons d’introduire dans la discussion 
va nous permettre de donner au probléme de la détermination des con- 
stantes de la courbe une autre présentation que celle que nous con- 
naissons déja. Il] faut que le procédé de calcul de ces constantes soit 
tel qu’il donne le méme résultat, lorsqu’il porte sur les mesures corre- 
spondant 4 un méme stade, et lorsqu’il est appliqué a deux ou a 
plusieurs stades successifs. Pour nous placer dans un cadre plus con- 
eret, nous adopterons un mode de transcription graphique des résultats 
qui rende sensiblement normal le nuage des points figuratifs de la 
population des mesures relatives 4 un stade donné. Nous savons que la 
résultat sera atteint fréquemment en utilisant les coordonnées log- 
arithmiques; parfois pourront étre employées tout aussi bien des 
coordonnées arithmétiques. La corrélation entre les deux variants 
etant forte par hypothése, le nuage sera trés allongé suivant une droite 
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sur Jaqueile s’échvlonneront aussi, d’aprés ce qui vient d’étre dit, les 
nuages figuratifs des autres stades. La définition des paramétres de 
cette droite doit étre telle que le résultat obtenu soit le méme, que 1’on 
ait utilisé un seul des nuages de points, plusieurs d’entre eux, ou leur 
ensemble. Un calcul trés simple montre que cette condition suffit 
exclure les droites de régression, que nous savions déja ne répondre 
que trés imparfaitement au probléme posé, et, au moins lorsque les 
coordonnées sont logarithmiques, l’axe principal qui, fournissant une 
solution symétrique en X et Y, aurait pu 4 la rigueuer étre retenu. I] est 
clair, en revanche, que si le rapport des variances de Y et X a méme 
valeur pour toutes les populations, ce que nous avons supposé, en ad- 
mettant que l’allométre est la méme pour deux stades consécutifs, le 
résultat cherché sera atteint en prenant pour pente de la droite le quo- 
tient des écarts types des variables comparées qui, si les nuages de points 
sont réellement alignés, est le méme pour chacun d’eux et pour le nuage 
global. Nous retrouverons ainsi, par une toute autre voie, le résulta‘ 
auquel nous étions arrivé en étudiant |’allométrie de taille. 
Ainsi se trouve établie, pour un Arthropede, la proposition an- 
noncée, mais avant de poursuivre, il importe de préciser que les con- 
clusions auxquelles nous sommes arrivés offrent un degré assez in 
de généralité. 


Nous avons démontré qu’a l’intérieur d’une étape de la croissance 
une méme loi f convient 4 la fois & 1’allométrie de taille pour chacune 
des stades et & l’allométrie de croissance pour l’ensemble de ces der- 
niers. Nous avons démontré accessoirement que deux étapes successives 
sont normalement séparées sur la courbe figurative du phénoméne par 
un point anguleux ou une discontinuité, marques du passage d’une loi 
de croissance 4 une autre. Mais nous n’avons pas démontré qu’il 
existait de telles étapes et ne pouvions le démontrer puisqu’il s’agit la 
d’une question de fait. L’observation montre d’ailleurs que s’il existe 
de nombreux cas ot une étape comporte cing ou six stades ou méme 
davantage, il en est d’autres, tels que les développements postpubéraux 
de certains Crabes, ot chaque stade doit étre considéré comme consti- 
tuant a lui seul une étape distincte. La courbe totale peut ainsi com- 
porter a son extrémité un ou plusieurs segments ne représentant qu’une 
allométrie de taille, alors que les précédents descrivent a la fois une 
allométrie de croissance et une allométrie de taille. C’est le cas pour 
Maia squinado, ot les courbes traduisant le développement des ap- 
pendices dans les deux sexes, ou celui de l’abdomen des femelles, com- 
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portent dans la partie actuellement connue, trois segments, dont le der- 
nier, séparé du précédent par une discontinuité, se rapporte exclusive- 
ment a l’allométrie de taille de 1’adulte. 

Plusieurs hypothéses supplémentaires ont été nécessaires pour que 
puisse étre précisée la forme de la relation f(z, y) =0. Nous avons 
supposé que nous disposions d’un mode de transcription de nos ré- 
sultats qui rendait sensiblement normal, 4 chaque stade, le nuage de nos 
points figuratifs. La démonstration précédente implique bien que, si 
ce résultat pent étre atteint pour |’un des stades, il l’est nécessairement 
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Allométrie de 1’abdomen chez les Maja squinados femelle, & deux stades 
eonsécutifs du développement = 1,28). 
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pour tous ceux qui appartiennent 4 la méme étape, puisque la relation 
qui unit les deux variables est unique et que devenue linéaire dans une 
de ses parties, elle le devient nécessairement dans les autres. Mais elle 
ne peut pas nous dire quel est le type d’anamorphose qui convient 4 
chaque cas. S’il est vrai que le résultat peut étre souvent atteint en 
utilisant des coordonnées logarithmiques, et si nous avons de bonnes 
raisons de croire qu’il en est trés généralement ainsi, nous connaissons 
quelques cas oU il en va autrement. Ce sont ceux qui sont connus sous 
le nom d’allométrie variable et ot la courbe logarithmique présente 
une convavité plus ou moins marquée ou méme un point d’inflexion. 
Ils s’observent surtout, semble-t-il jusqu’a présent, sur les trongons de 
la courbe totale qui se rapportent aux adultes et représentent essen- 
tiellement des allométries de taille. Mais il n’est pas exclu que l’on 
puisse trouver quelque jour des allométries de croissance du méme 
type. 

Enfin nous rappellerons que le choix de la ‘‘meilleure courbe’’ repré- 
sentative d’un ensemble de mesures comporte toujours une large part 
d’arbitraire et que la solution 4 laquelle nous sommes arrivés, pour 
satisfaisante qu’elle puisse étre, n’est pas la seule que 1’on soit en droit 
de retenir. Nous préciserons également que cette solution est essen- 
tiellement théorique, qu’elle n’est valable que lorsque sont remplies au 
préalable toute une serie de conditions et qu’au surplus elle ne prétend, 
en aucune maniére, 4 étre une méthode réguliére de calcul des cons- 
stantes d’une courbe d’interpolation. , 


Les remarques et les réserves qui viennent d’étre faites, pour utiles 
qu’elles puissent étre, ne doivent cependant pas étre considérées comme 
diminuant la portée de notre démonstration. Elles en précisent seule- 
ment la signification et rappellent, s’il en était besoin, qu’un schéma 
théorique ne doit pas étre appliqué 4 un cas concret sans d’indispen- 
sables précautions. I] nous reste 4 montrer rapidement que les ré- 
sultats précédents s’étendent aux organismes dont la croissance ne se 
fait pas par mues. Si nous avons choisi, pour notre démonstration, le 
cas des Arthropodes, ¢c’est parce que, dans cet embranchement, les 
divers stades définissent clairement une succession de populations sur 
chacune desquelles il est possible d’étudier les phénoménes d’allométrie 
de taille. Mais il est clair qu’une telle étude peut également étre faite 
sur des animaux 4 croissance continue, si l’on posséde un critére per- 
mettant de définir sans ambiguité les individus comparables, critére 
qui évidemment pourra la cas échéant étre l’age. Les raisonnements 
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précédents restent valables, sous la réserve habituelle que, dans la pe- 
riode étudiée n’intervienne pas de phénoméne physiologique capable de 
donner 4 deux animaux de méme taille des structures dissemblables. 
Pour le reste, la différence intéressante entre la cas général et celui des 
Arthropodes porte sur la définition des stades successifs qui, au lieu 
d’étre imposée par la nature méme du développement, peut devenir 
assez largement arbitraire. On peut, en principe, multiplier ces stades 
autant qu’on le désire, sans que rien soit changé 4 nos raisonnements, 
ce qui permet d’étendre la démonstration du cas des animaux 4 dé- 
veloppement discontinu 4 celui des animaux 4 croissance continue. 

Quelques lignes nous suffiront pour étendre 4 la comparaison bio- 
chimique des organismes les procédés employés dans la comparaison 
morphologique. Rien, dans les formules que nous avons utilisées, 
n’implique que z ou y doivent mesurer plutét la longueur d’un appen- 
dice que le poids du cerveau ou la quantité de calcium renfermée dans 
un organisme. Toutes les interprétations valables dans un cas le sont 
aussi dans l’autre et la notion d’allométrie biochimique, qui a été vive- 
ment critiquée, mais qui a aussi été employée avee un grand succés, 
sans jamais pouvoir étre expliquée de fagon satisfaisante, apparait 
maintenam comme pleinement justifiée. 


Les nombreuses recherches de morphologie quantitative de ces 
vingt derniéres années ont permis de découvrir un grande nombre de 
faits nouveaux et de donner une forme précise 4 des notions dont les 
travaux antérieurs avaient pu faire pressentir l’importance. Ces re- 
cherches n’ont mis le plus souvent en oeuvre que des méthodes trés 
élémentaires et un appareil mathématique trés réduit. 

La simplicité de ce mode d’étude a beaucoup facilité sa diffusion ; il 
a permis d’acquérir en assez peu de temps des renseignements précieux 
sur la croissance d’animaux trés divers et d’établir des comparaisons 
précises entre adultes de méme espéce ou d’espéces voisines. II est 
certain que des recherches de ce genre conservent a l’heure actuelle 
une bonne part de leur utilité et que, longtemps encore, elles appor- 
teront, dans certains domaines, des contributions trés appreciables 4 
nos connaissances. 

Mais, dans un nombre croissant de circonstances, il apparait de plus 
en plus clairement que des conclusions légitimes ne pourraient étre 
tirées de telles études qu’a la condition de disposer de critéres plus 
sirs que ceux que permet le simple ajustement graphique. Aussi 
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a-t-on été conduit 4 joindre & la technique rapide et facile, mais seule- 
ment approximative, de l’allométrie, les techniques plus laborieuses et 
plus savantes, mais aussi plus rigoureuses, de la biométrie classique. De 
plus en plus fréquemment est apparue, au cété des constantes de la 
formule d’allométrie, l’indication de leurs écarts type. Ce faisant, on 
se contente d’ailleurs de juxtaposer l’emploi de deux procédés, 1’un 
d’origine biologique, 1’autre d’origine mathématique, sans lesefondre 
en une technique cohérente. 

La notion méme d’allométrie, malgré toutes les tentatives qui ont 
éte faites pour en pénétrer la signification biologique, n’est encore 
aujourd’hui guére plus qu’une simple donnée de 1’expérience, dont le 
fondement rationnel reste mal assuré, et que beaucoup de ses partisans 
considérent toujours comme essentiellement empirique. La diversité 
méme des circonstances dans lesquelles il parait permis d’en faire usage 
en rend la compréhension plus obscure et la contradiction s’aggrave 
entre l’importance grandissante des résultats obtenus et la précarité 
de plus en plus apparente de leur fondement. Pour surmonter cette 
contradiction, il est devenu nécessaire de repenser la notion méme 
d’allométrie. 

C’est ce que nous avons essayé de faire dans les pages qui précédent. 
Eliminant les difficultés accessoires, négligeant provisoirement tous les 
phénoménes qui, plus ou moins légitimement, apparaissent en désaccord 
avec la simplicité de la relation fondamentale, nous avons essayé de 
comprendre la signification véritable de la relation d’allométrie. Il 
m’apparait que nous y avons & peu prés réussi et qu’elle n’est rien 
d’autre qu’un aspect de la variabilité essentielle de tous les étres 
vivants. 


DISCUSSION 


Lois M. Zucker. I should like to discuss some data we have which 
fit in rather well with a question Dr. Teissier has raised. We have data 
on the relations between the femur ash and body weight in growing rats, 
selected in five age groups. Each age group provides a scatter diagram 
relating log-ash to log-weight and these scatter diagrams overlap as in 
Dr. Teissier’s example. Now, does the same line or the same slope 
apply both within ages and between ages? In other words is the rela- 
tion independent of age? The choice of line to use in fitting the data 
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becomes rather critical. Neither variable is independent and there 
is error in both variables. Therefore the preferred line may fall be- 
tween the two regression lines. If the deviations parallel to log-ash are 
minimized the five slopes do not differ significantly from a common 
trend which, in turn, is very much flatter than the trend between ages; 
if the deviations parallel to log-weight are minimized, the five slopes 
have a common trend whick does not differ significantly from the trend 
between ages. 

The line suggested by Dr. Teissier lies somewhere between. But is 
it the best line? It seems to us that it has a flaw. Consider the ideal 
line expressing the relationship. Suppose there are only, say, ten pos- . 
sible type positions along this line instead of a continuous set. Neither 
variable is independent. Actual points will range themselves around 
their type position in elliptical distributions which, accordingly to the 
error distribution, may be nearly circular (each variable with equal 
error) or tall and narrow (most of the error with the Y coordinate) or 
low and wide (most of the error with X). In the last case the regres- 
sion of X on Y is a good estimate of the true line while the regression of 
Y on X, and to a leaser extent lines like the one suggested by Dr. Teissier 
which he approximated midway between the two regression lines, are 
too flat in slope. In the case with most of the error with the Y coordi- 
nate, the regression of Y on X isa good estimate, where the regression of 
X on Y, and to a lesser extent Dr. Teissier’s line, are too steep. It is 
only if the error is approximately the same for both variables that a line 
between the two regression lines, such as that suggested by Dr. Teissier 
is a good estimate. It does not seem possible to accept a single solution 
independent of the error distributions for all these situations, like Dr. 
Teissier’s line or the one suggested by Dr. Wald in a recent publication. 

In response to Dr. Teissier’s suggestion that his line is recommended 
only when neither variable is independent, I should like to state that 
a scatter diagram like that obtained with one variable independent can 
also arise for neither variable independent, if the error is mostly associ- 
ated with one variable. It is usually impossible to estimate objectively 
the relative error associated with the two variables in cases like this, 
because the error is largely biological and not error of measurement. 

We believe, for instance, that in the relation between femur ash and 
body weight much more error is associated with body weight than with 
ash. It is very difficult to affect the ash by experimental means—water . 
deprivation, starvation, disease, or nutritional deficiencies other than 
of calcium, phosphorus or vitamin D, while the body weight, largely 
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affected by the weight of water and soft tissues, is very sensitive to ex- 
perimental procedure. Should it not, therefore, be much more vari- 
able? For this situation we should prefer a line very close to the re- 
gression of log-body weight on log-femur ash, not a line halfway to- 
wards the other regression line. However, our belief as to the error 
distribution cannot be directly verified objectively. 

(For data see Zucker and Zucker, American Journal of Physiology, 
146, 585, 593, 1946 ; for further discussion see Zucker, Human Biology, 
19 : 232, 1947.) | 


Harold Hotelling. From a purely descriptive standpoint, data on 
a pair of variates such as measures of a bone cannot possibly be sum- 
marized adequately by any one regression line. Even if we make the 
simplifying assumption of a bivariate normal distribution, there are 
five parameters, and these cannot be summarized by the two coefficients 
in the equation of a straight line. 

When we pass over to the relation of bone dimensions or chemical 
composition to vitamin D in the diet we are on different ground. If 
the object is simply to predict bone ash percentage as a function of 
vitamin D, then asimple regression equation for bone ash percentage on 
vitamin D quantity is in order, and this equation will disregard all 
measures of the bone excepting ash content. If, however, there is a 
question whether vitamin D has some effect on the bone, without ad- 
vance knowledge of the particular kind of effect, and if say eight meas- 
ures on each bone are taken to represent it fully, then a more general 
method should be used. This is arithmetically identical with the tech- 
nique of fitting a multiple regression equation and testing the sig- 
nificance of R. 

The question at issue may be the still broader one of finding some 
combination of vitamins that will affect in some way, not initially 
known, the shape or composition of a bone. Then the appropriate 
method is to determine together the most predictable function of the 
bone dimensions and contents, and the best function of the various 
vitamin doses for predicting this. Techniques for doing this, and for 
testing the significance of the results, have been developed (Biometrika, 
28 : 321-377, 1936). This kind of approach will, I think, supersede, for 
this kind of purpose, the summarization of the relations between meas- 
ures on a bone by any regression line or its n-dimensional generaliza- 
tions. 
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N. Rashevsky. The statistical relations discovered by Professor 
Teissier constitute important new knowledge. This kind of knowledge 
concerns how phenomena happen. If we wish the answer to the ques- 
tion of why they happen, we must extend our approach by introducing 
hypotheses and postulates. Professor Teissier is well aware of this 
himself, since he mentioned that he had a physiological explanation of 
some of these phenomena. 

To illustrate what the theoretical method can do with the problem of 
growth, I wish to mention some studies in the Section of Mathematical 
Biophysics at the University of Chicago, relating to the problem of can- 
cer. 

Multicellular organisms grow up to a certain limit, whereas in tissue 
cultures cellular multiplication may go on practically indefinitely. 
This suggests that some inhibitory factor is produced by each cell, so 
that as the total number of cells of the organism increases and each cell 
becomes inhibited by a greater and greater number of cells, growth 
eventually stops. The mathematical formulation of such an assump- 
tion leads to growth curves reminding us of those described by Professor 
Teissier. As has recently been shown by Kesselman, some aspects of 
allometric growth can be described similarly. The theory implies that 
after the adult stage is reached, the total amount of the inhibitory factor 
decreases due to the repair of the natural wear of the organism and of 
accidental wounds. The theory leads to an expression for the decrease 
of velocity of wound-healing with age. Although the data on this ques- 
tion are very meager, whatever is available is in agreement with the 
theory. As the organism grows older and the total inhibitory effect 
decreases, small accidental fluctuations will be sufficient to produce 
accidental growth. Hence such accidental growth would be expected 
to be more frequent in older age. Actually, it is known that the inci- 
dence of cancer does increase with age. 

The theory also enables us to derive the incidence curve as a function 
of age and the theoretical curve is found in good agreement with the 
observed one. From one of the parameters of the incidence curve we 
can compute the natural human longevity, longevity as it would be if 
unaffected by illness or accident. The theoretically computed value 
turns out to be of the order of 80 to 100 years. It is certainly remark- 
able that from a single theory we may compute such different things as 
decay in wound-healing ability, incidence of cancer and human longev- 
ity. But the most important progress in science consists principally 
in finding relations between phenomena which at first sight appear to 
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be unrelated. In cancer research the purely statistical method is of 
extreme importance to give us the data. But for a complete success it 
must be supplemented by a broader theoretical approach. 


Georges Teissier: J’ai suivi avec intérét les remarques prégentées 
4 la suite de mon exposé, mais si je crois avoir compris le sens général 
de ces interventions, je ne suis pas assez assuré d’en avoir suffisamment 
saisi les nuances pour pouvoir les discuter en détail. 


A. Mme. Lois ZucKEr je repondrai seulement qu’il ne me semble 
pas que nous parlions exactement des mémes choses. Les hypothéses 
dans lesquelles elle se place, pour les deux courbes qu’elle étudie, diffé- 
rent trop des miennes pour que des conclusions valables dans le cas 
dont je me suis occupé puissent s’appliquer sans modification aux 
schémas qu’elle a envisagés. 


Au Dr. Hore.uine, je répondrai que je sais trés bien que la ligne 
d’allométrie, non plus qu’aucune ligne de régression d’aucune espéce 
ne peut fournir 4 elle seule la description compléte d’un ensemble de 
mesures. Je n’entends pas substituter 4 la description statistique 
classique des faits un procédé plus simple et aussi efficace, mais bien 
fournir de 1’essentiel de ces faits, un schéma clair et interprétable en 
termes biologiques. L’immense litterature consacrée depuis plus de 
vingt ans & la croissance relative établit sans conteste que la relation 
d’allométrie remplit parfaitement ce but. Je me proposais au 
jourd’hui, non pas d’établir sa validité pratique, qui est certaine, mais 
seulement de dégager sa double signification, biologique et statistique, 
résultat auquel je crois bien étre parvenu. 


A M. RasHevsxy je dirai que je suis moi-méme responsable d’un 
schéma physiologique qui permet, a ce que je crois, de justifier de fagon 
trés suffisante l’application de la relation d’allométrie 4 un organisme 
en croissance, et qui explique en méme temps |’existence des stades 
critiques et des étapes du développement. Mais ce schéma ne vaut ni 
pour l’allométrie de taille, ni pour 1’allométrie biochimique. 
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Comme il serait déraisonnable de supposer que 1’identité des rela- 
tions qui intervienrent dans les trois cas puisse résulter d’une simple 
. coincidence et que la justification d’une relation applicable en des cir- 
constances aussi dissemblables puisse reposer sur des considérations 
mécanistes, il faut chercher dans une autre voie. Je crois avoir réussi a 
montrer qu’une interprétation statistique répond entiérement a la 
question posée. Cela ne signifie naturellement pas qu’il faille renoncer 
4 connaitre les mécanismes physiologiques et physicochimiques qui in- 
terviennent dans chaque cas. Tout en étant décidé, pour ma part, a 
continuer l’étude que j’ai entreprise de cet aspect du probléme, je 
tiens 4 déclarer qu’A mon avis, et tenu compte d’innombrables tenta- 
tives dont celle de Robertson reste le prototype, il est dangereux 
de vouloir fonder une théorie de la croissance sur un schéma physico- 
chimique nécessairement trop simple. 
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THE GENERAL THEORY OF PRIME-POWER 
LATTICE DESIGNS* 


I. INTRODUCTION AND DESIGNS FOR p* VARIETIES IN 
BLOCKS OF p PLOTS 


Oscar KEMPTHORNE AND WALTER T. F'EDERERT 


I. INTRODUCTION 


Extensive use is now made of the lattice designs originated by Yates 
[5, 6] for testing large numbers of varieties. The general principle in 
the structure of these designs is that the varieties are regarded as aris- 
ing from the combination of factors at the same or different levels. 
Some of the effects and interactions between pseudo-factors are con- 
founded with blocks or with some other restrictions in each replicate 
of the experiment. 

A number of replicates must be used, confounding, in general, dif- 
ferent sets of effects and interactions. Such confounding results in a 
gain in information on unconfounded effects and interactions over what 
would be obtained with randomized complete blocks, since they are 
based on comparisons within smaller blocks of plots. This gain may be 
offset to some extent by a loss in information on effects and interactions 
which are confounded among incomplete blocks. 

Yates [8,9] and Cochran [1, 2] have described the theory by which 
information contained in block comparisons can be utilized for the 
various designs which had been devised. The purpose of the present 
series of papers is to give a systematic description of lattice trials with 
any number of replicates for a number of varieties which is a power of 
a prime. The present paper contains the basic factorial theory and a 
description of designs for p” varieties in blocks of p plots. It will be 
followed by papers dealing with other lattice designs and numerical 
examples. 


* Contribution of the Statistical Section of the Iowa Agricultural Experiment Sta- 
tion in cooperation with the Bureau of Agricultural Economics, United States Depart- 
ment of Agriculture. Journal paper no. J 1553. Project 890. 

t Associate Professor, Statistical Laboratory, Iowa State College, and Associate 
Agricultural Statistician, Bureau of Agricultural Economics, ommanaenting with the Iowa 
Agricultural Experiment Station, respectively. 


54 


| 
| 
t 
| 
aa 
| 
VE 
i 
{ 
| 
4 
vin 
4 
= 


PRIME POWER LATTICE DESIGNS 


Il. BRIEF DESCRIPTION OF THE p" FACTORIAL SYSTEM 

(a) Effects and interactions. A description [3] has been given of 
the p* factorial system (where p is a prime) with the use of geometrical 
terminology and it is sufficient here to recapitulate the main results. 
The p" treatment combinations may be represented by an n-dimensional 
lattice, each side of which contains p points: thus vsing coordinates 
Zi, Za, . . . » L_ to represent the levels of each factor in a treatment com- 
bination, the control-treatment combination will be represented by 


the treatment consisting of the first factor at unit level and all the other 
factors at zero level by 


.=2,=0; 


and the treatment combination in which all factors are at level (p—1) 


by the point 


Using the ordinary definition of effects and interactions as given in 
Yates’ fundamental work, ‘‘The Design and Analysis of Factorial Ex- 
periments’’ [7], the main effect of factor 1 is given by the contrast 
between the yields of treatment combinations represented by 


z,=0, z,=1, 2, =2, 


and it has (p—1) degrees of freedom. There are in all n such compari- 
sons for main effects which may be represented symbolically by 21, 22, 
Zs, ..., Zn, each having (p—1) degrees of freedom. 

The interaction of two factors represented by x, and 2, has in all 
(p -1)? degrees of freedom and these may be split into (p-—1) sets of 
(p—1) degrees of freedom represented as follows: 
the contrasts between those treatment combinations for which 


21 +2,=0, 1, 2, p-1 (mod p), 
the contrasts between those treatment combinations for which 
2, + 22, =0, 1, 2, p-1 (mod P), 


and so on to the contrasts between those treatment combinations for 
which 

+ (p—1) 1, 2,..., p—1 (mod p). 
The formal relationship of the above definitions of interactions with 


: 
: 
| 
= 
| 


56 BIOMETRICS MARCH 1948 


those given by Yates [7] in the case of the 3" system is easily seen. 
The four degrees of freedom for the interaction of two factors a and b 
each at three levels are, in Yates’ terminology, partitioned into two 
pairs of degrees of freedom, AB(I) and AB(J), each pair being ob- 
tained from the comparison of three totals. By reference to Yates’ 
definitions [7], the contrast AB(J) is the same as that denoted above 
as the contrast between the three totals given by 


=0, 1, 2 (mod 3), 
and that AB(J) is the same as that denoted by 
2, + 2x, = 0, 1, 2 (mod 3). 


The three-factor interaction of z,, z, and z; has in all (p—1)* 
degrees of freedom and may be split into (p—1)? comparisons of p 
totals, each comparison having p—1 degrees of freedom. In our geo- 
metrical terminology they are given by 


+ + = 0, 1, 2, p-1 (mod p) 


where a, and a; take on all values from 1to (p-—1). Again the formal 
identity with Yates’ definitions [7] is easily seen for the case of three 
factors each at three levels : in this case, thé degrees of freedom for the 
three-factor interaction may be divided into four sets of two degrees of 
freedom, represented by 


1, 2 (mod 3), 
2_,+2z, =0, 1, 2 (mod 3), 
+22,+2,=0, 1, 2 (mod 3), 
+ 22, + 22; =0, 1, 2 (mod 3), 


and these are the comparisons which were denoted by Yates as Z, Y, X, 
W respectively [7]. 

The extensions of the above are obvious and need not be given in 
detail. If the factors are called a, b,c, d etc. then with a simplification 
of the notation the effects and interactions may be represented as 
follows : 
main effects A, B, C, D, ete. 
two-factor interactions 


AB, AB?,... , AB?-}, 
BC, BC?,.. . , BC?-, ete., 
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three-factor interactions 


ABC, ABC?, AB*C, AB*C?, ... , AB*-* C?-}, 
ABD, ABD*, AB*D, AB*D*, . . . , AB*-* D®-', 
etc. 


It should be noted that in order to have a unique enumeration, the 
power of the first symbol should be unity. This may be done by noting 
for example that A*BC gives the same contrasts as AB*C?. 

In the 2" factorial system, the concept of generalized interaction 
[7] is of great use: If, for example, interactions ABCD and BCDE are 
confounded then so is their generalized interaction AE, obtained by 
multiplying the symbols together with the rule that 


A? = .=1. 


Similarly, the concept of generalized interaction is necessary for the 
discussion of the p* factorial system. It may easily be shown by use 
of the geometrical method described above that if effects or interactions 
represented by X and Y are confounded, then the effects or interactions 
represented by XY, XY*,. .., XY*-* must also be confounded. By 
the generalized interaction of two interactions A“ B"C™ and 
. . . is then meant the (p—1) interactions 


ASH 
Bhi + 2% 


where (i) the powers are all to be reduced modulo p (that is divided 
by p and the remainder substituted), or as is the same thing where use 
is to be made of the relation A’ = B?=C* . . . =1, and (ii) where the 
power of the first letter or any symbol is forced to be unity, by taking 
when necessary the power of the symbol which will make this the case. 
To give a concrete example, in the 3° system with factors a, b, and c, 
two of the four three-factor interaction pairs of degrees of freedom are 
ABC and ABC’: their generalized interaction consists of 


ABC - ABC? = A*B?C* = A*B? = (A*B?)? = A*B* = AB, 
and 
ABC - A?B*C* = A®B°C' = C2 =C. 


If then a 3° factorial experiment is arranged in blocks of three plots 
so as to confound ABC and ABC’, AB and C will also be confounded. 
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(b) The yields of treatment combinations in terms of effects and 


‘interactions. In the case of the 2" factorial experiment, it is well 


known that the yield of a treatment combination may be expressed in 
terms of effects and interactions. Thus if three factors, a, b, c, are 
tested, the yield of the treatment combination a, b; cx is given by the 
expression : 


mean + (—1)*-? 4 + (—1)’-? B+ (—1)*"? C+ (-1)**/ AB+ 
(—1)*+® AC + (—1)!** BC + (—1)**#**-1 ABC], 
where the effects and interactions have been reduced to a single-plot 
basis. 
A is here defined as the difference between the mean yield of plots 


receiving factor a at the unit level and the mean yield of plots receiving 
the factor at the zero level; or in algebraic terms 


A =}(a-1)(b+1)(ce+1), 


where the algebraic expression is to be expanded and the correspond- 
ing yields substituted for the treatment combinations; the interaction 
AB is likewise defined as 


AB = 3(a-1)(b-1)(c+1), 
and so on. 


In the general case with factors a, b,c, . . . each at p levels, effects 


or interactions cannot be represented by a single difference, since there 
are in all (p—1) differences. In this paper, the mean yield of the plots 
for which factor a was at level i is denoted by (A),. In the case of 
interactions we denote by (A*B®C”) at + Bj + yk, the mean yield of those 
plots for which 

aX, + BX, + = at + Bj + yk (mod p). 


which is therefore the mean yield of the plots that enter into the partic- 
ular component of the interaction ABC, specified by the plot a,bjcx. 
Thus, considering the plot in the 3* system which receives the treatment 
combination a,b,c, and the interaction ABC?, 


at, + Bi, = 1(1) +1(2) +2(1) =5=2 (mod 3), 


and the quantity (A°*B°C”) ai + 8j + yk becomes (ABC*), and is defined 
to represent the mean yield of those plots for which 


%1+2%,+22,=2 (mod 8). 
With this notation, the yield of the treatment combination ajb;cx . . . 
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+ (AB) 4,5 .. + (ABP?) 
+ (AC) (AC?) (AC?-*) $+ (p-1)k 
+ (ABC) 4,5,%+ + (ABP? 4, 


where the suffices are all reduced modulo p. 
That this gives the same as the above formula for the 2" system may 
be seen by noting that: 


A = (A),—-— (A), =2[(A)1— mean] 
and so on. 


Il. TRUE FACTORIAL DESIGNS AND QUASI-FACTORIAL DESIGNS 


The above formulation of the factorial scheme is particularly useful 
in the consideration of true factorial designs. In the particular case 
when p = 2, the formulation is identical with that which has always been 
used to devise systems of confounding and fractional replication. In 
the case when p = 3, the present formulation may be used with advan- 
tage to simplify the discussion of confounding and fractional replica- 
tion for the system. 

In the true factorial system, when n individual factors are tested, 
use is made of the fact that main effects and interactions between a 
small number of factors are of importance, but that interactions be- 
tween many factors are likely to have negligible values in comparison 
with the experimental error. It is therefore possible to arrange the 
experiment in blocks, confounding between blocks interactions which 
the experimenter considers unimportant, with a resultant increase in 
precision. 

In the quasi-factorial design, on the other hand, the effects and inter- 
actions are purely formal and all are of equal importance, as the table 
of variety means which is finally required can only be obtained with 
knowledge of all effects and interactions, no matter what order. It is 
_‘ Mecessary to adopt some formal confounding within each replicate in 
_' oder to reduce the experimental error by reducing the size of block. 

Within any one replicate, some of the effects or interactions are 
obtained by comparison of block (or row or column) totals and the re- 
mainder by comparisons within blocks. The comparisons of block 
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totals is subject to an error denoted by the inter-block error, and com- 
parisons within blocks are subject to the intra-block error. In general, 
the inter-block error is considerably greater than the intra-block error 
[4]. In order that estimates of reasonable precision may be made for 
all the effects and interactions it is necessary that different effects and 
interactions be confounded in the different replicates. 


IV. DESIGNS FOR p" VARIETIES IN BLOCKS OF Pp PLOTS 

An outline of the possible designs, considering any number of repli- 
cates is now presented. 

(a) m=2. Ifthe factors are denoted by a and b, then arrangements 
in blocks of p plots exist which confound one of the following effects or 
interactions: A, B, AB, AB?,..., AB*’-*. The case in which A is 
confounded in one replicate and B in another replicate is called the 
two-dimensional lattice with two sets, or the simple lattice. A design 
in which A is confounded in one replicate, B in another, and one of 
AB, AB?*, AB*,... , AB*- ina third is called the triple lattice. It 
is possible to have p+1 replicates in all, when each effect and inter- 
action is confounded in one of the p+1 replicates and unconfounded 
in the other p replicates. 

(b) m=3. Let the factors be a,b,c. It is possible to use blocks of 
p plots, with p? blocks randomized within each replicate, the effects and 
interactions which are confounded with blocks being chosen from A, B, 
©, AS, AO, « OO, ..., 
ABC, ABC?,..., AB*C,..., AB*-?C*-*. The group of 
interactions confounded in any one replicate will consist of two effects 
or interactions with their generalized-interaction as defined above. If 
for example the effects A, B are chosen to be confounded in a replicate, 
then so must the interactions AB, AB’, ... , AB*-*. It is clear that 
several choices are possible and later in this paper particular examples 
are discussed. If, as is generally the case, intra-block information is 
required on all effects, at least three replicates must be used. 

(c) General n. Let the factors be denoted by a,b,c,.... Then 
a replicate in p"-* blocks of p plots could be obtained by confounding, 
for example, all the effects and interactions of (p—1) of the factors. 
The principles on which designs are based will become clear from the 
consideration of particular cases. In general, at least n replicates must 
be used. This paper does not deal with designs for a number of vari- 
eties which is the power of a non-prime, but it will be clear from the 
discussion that n-dimensional lattices with n replicates of k*-* blocks 
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of k plots for a number of varieties k", where k is not a prime, may be 
constructed by confounding each pseudo-factor in all but one of the 
m replicates. The number of different replicates possible depends on 
the properties of kx k Latin Squares. Some of the properties are 
discussed by Finney [10]. 

(d) p® system or 3-dimensional lattice in blocks of p. (i) Design. 
Let the factors be denoted by a, b, c, each having 3 levels for purposes 
of illustration. Then it is necessary in each replicate to confound 
8[= (p?—1)] degrees of freedom. The possible schemes of confounding 
are given in Table 1. 

The order in which the replicates are given in Table 1 is logical 
and arises as follows. Consider first effect A to be confounded; if B 
is also confounded, so then are AB and AB’, and if either AB or AB? is 
confounded then so is the other and B; next, the group A, C, AC, AC? 
is found : the next possible interaction to be confounded along with A is 
BC and this results in the confounding of ABC and AB?C?: finally 


TABLE 1 
PossIBLE CONFOUNDING SYSTEMS FOR 33 SYSTEM IN BLOCKS oFr 3 


Effect or Interaction 


System | A B AB AB? | C | AC | AC? | BC | BCt | ABC | ABC? | | AB°C: 
1 x x x x 
2 x x x x 
3 x x x x 
4 x x x 
5 x x 
6 x x x 
7 x x 
8 x x 
9 x x x 
10 x x 
11 x x x x 
12 x x x 
13 x x x x 
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if BC? is confounded along with A then so are ABC? and AB?C: all the 
main effects and interactions have then been confounded with A in one 
of the replicates. Consider next the effect B to be confounded: it has 
already been considered with effect A and interactions AB and AB? so 
we may proceed from C across the table as with A previously. In this 
way the 13 (in general p? + p+1 for a p* system) possible systems of 
confounding of the 3° system in blocks of three have been generated. 
The constitution of the block for a particular system of confounding 
may be determined from consideration of the definition of the effects 
and interactions. The constitution of the blocks on the case when A, 
B, AB and AB? are confounded is simply obtained. A less obvious case 
is considered here, say when AB*, AC BC and ABC? are confounded. 
Let the varieties be denoted by a,b,c, where each i, j, k run from 0 to 2. 
Then the combinations which enter into groups of three blocks are as 
follows : (the suffices t, j, k are used to denote the particular varieties. ) 


$+2j = 0 (mod 3) : 000, 001, 002, 110, 111, 112, 220, 221, 222 
AB?* {4 +2j = 1 (mod 3) : 020, 021, 022, 100, 101, 102, 210, 211, 212 
$+2j = 2 (mod 3) : 010, 011, 012, 120, 121, 122, 200, 201, 202 


++k =0 (mod 3) : 000, 010, 020, 102, 112, 122, 201, 211, 221 
AC{i+k =1 (mod 3) : 001, 011, 021, 100, 110, 120, 202, 212, 222 
t+k =2 (mod 3) : 002, 012, 022, 101, 111, 121, 200, 210, 220 


It is easy to see that the nine blocks are: 


I VI VIE VII IX 
000 «4001 002 «(020 «#4021 «#4022 #+#O10 012 
112 «#110 #111 #4102 #4100 101 £122 120 121 
221 222 220 211 212 210 201 202 # 200 


The contrast represented by AB? for example is then given by the con- 
trast of blocks I+ II + III, blocks IV +V+VI, and blocks VII + VIII 
+ XI. 

In planning variety trials of the type considered in this paper for a 
particular number of replicates, it is important to choose systems of 
confounding (as given, for example, in Table 1) so that the total amount 
of confounding is spread as far as possible over all the possible effects 
and interactions. The information available on any particular effect 
or interaction will be of two types (a) intra-block, that is, based on 
within-block comparisons and (b) inter-block, that is, based on com- 
parisons among block totals. In general, information of type (a) will 
be more accurate since block variability, measured by differences among 
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blocks after taking account of variety differences, will be large com- 
pared with the within-block variability [4]. The experiment is there- 
fore designed to yield as much information as possible of type (a). 

If three replicates are to be used, a number of choices of equal value 
are available which, by permutation of symbols, are equivalent to the 
choice of replicates consisting of numbers 1, 2, and 5: that is, with 
A, B, AB and AB* confounded in one replicate; A, C, AC and AC? con- 
founded in a second replicate; and B, C, BC and BC? confounded in a 
third replicate. If four replicates are to be used, the best choice con- 
sists of the following: 


number 1 confounding A, B, AB, AB? 


number 2 7 A, C, AC, AC? 
number 5 = B, C, BC, BC? 
number 9 AB, AC, BC?,  AB*C? 


Replicate number 10, 12, or 13 may be used in place of number 9. 
With this choice, three main effects and three interactions are con- 
founded in two out of the four replicates, and four interactions are 
confounded in one of the four replicates. The remaining three inter- 
actions are completely unconfounded. 

In the same way the best choice for any number of replicates up to 
13 may be considered. The systems of confounding are rearranged in 
Table 2 so that the best design (or the more nearly best) for r replicates 
is given by the first r rows of the table. The table was constructed on 
the principle that all the effects and interactions should be confounded 
in as equal a number of replicates as possible: on investigation it was 
found that there is no unique order but the order om is nearly the 
best under all circumstances. 

(ii) Analysis. A method for the k* system with 3 replicates con- 
sisting of numbers 1, (or X), 2, (or Y), and 5, (or Z) or multiples of 
these three replicates has been described by Yates [8]. It is not the 
purpose here to give a detailed computational procedure for all the 
possible cases, but merely to indicate the structure of the analysis and 
of the procedure for estimating adjusted varietal means. The analysis 
is given in terms of the concepts used in this paper. The confounding 
was as follows: 

replicate X: A, B, AB, AB’, 
Y: A, C, AC, AC’, 
Z: B, C, BC, BC. 
The structure of the analysis of variance is given in Table 3: 
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TABLE 2 


BIOMETRICS 


ScHEMES OF CONFOUNDING IN ORDER OF VALUE 
Effect or Interaction 


MARCH 1948 


A AB | C | AC BC | BC? | ABC] ABC: AB*C | AB*C? 
1 x x x 
2 x x x 
5 x x 
9 x x x x 
13 x x x 
12 x x x 
10 x x x 
11 x x x 
x 
3 x x x x 
e x x x x 
7 x x x 
« x x x x 
TABLE 3 
STRUCTURE OF ANALYSIS OF VARIANCE FOR 3 REPLICATES 
oF p3 SYSTEM IN BLOCKS OF p PLOTS 
DF Expectation inf Mean Square 
General p=3 General 
Replicates 2 2 
Blocks component (a) 8(p-1) 6 os? + pov? 
component (b) 8(p-1) + 1/3por* 
component (c) 3(p-1)? 12 + 2/3pon* 
Varieties (ignoring blocks) pt-1 26 
Error (intra-block) 2p* — 3p?+1 28 oi? 
Total 3p?-1 80 
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The degrees of freedom and sum of squares for blocks component (a) is 
obtained by noting that the blocks of replicates X and Y may both be 
grouped to estimate the effect A; the interaction of these two replicates 
and the A effect contributes two degrees of freedom: comparison of 
effect B in replicates X and Z, and of effect C in replicates Y and Z each 
contribute two degrees of freedom, making a total of 6. The blocks 
component (b) is obtained by comparing the mean A effect in replicates 
X and Y against the unconfounded A effect in replicate Z, and corres- 
pondingly for effects B and C, giving 6 degrees of freedom. The 
blocks component (c) is obtained by comparing each interaction which 
is confounded in one replicate only with the unconfounded interaction 
estimated from the other two replicates; each comparison yields two 
degrees of freedom and the six interactions altogether give a total of 12 
degrees of freedom. The sum of squares for varieties ignoring blocks 
is computed directly from the variety totals, and the sum of squares for 
replicates in the usual way. The error term is obtained by subtraction. 

The expectations of the mean squares of interest as stated by Yates 
[8], are given in the last column of Table 3 in terms of the within-block 
variance o,? and the additional between-block variance o»?. 

The evaluation of effects and interactions for this case proceeds 
exactly as described by Yates [8], but since only effects and two factor 
interactions are confounded it is convenient to state the method in a 
more general form. Consider the effect A for example: it is con- 
founded in replicates X and Y but is unconfounded in replicate Z. 
Estimates may be made for each of the replicates to make up a table of 
the form : 


(A)o (A), (A). 
Rep. X (A)ox (A)ix (A) ex 
Y (A) or (A)iy (A)ey 

Z (A) oz (A).z (A)az 


For the comparisons of (A) , (A), and (A). in replicates X and Y, 
which are inter-block comparisons, each quantity is subject in the gen- 


eral case to variance | Poy* + a” {| while for replicate Z the compari- 


sons are intra-block, each quantity having a variance o;7/,’. , If two 
quantities a and 8 are independent estimates of a parameter with vari- 
ances o” and o””, respectively, the linear function of a and 8 which is 
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the best estimate of the parameter in the sense of having the least pos- 
sible variance is 


wo+w’B 
where w = 5 and w’ = Gz» and this estimate will have a variance of 
. In the present case putting equal to w 
w+w’ 4 + 
w’(A)ox+w’(A)oy +w(A) oz 


the best estimate of (A), is and likewise 


2w’ +w 
for (A), and (A)... The variance of each of these estimates for com- 


parisons among themselves is In the simple case when only 


2w’+w- 
main effects and two-factor interactions are confounded, it is simpler 
to put the best estimates in the following form, exemplifying that for 
(A)o, 


indicating that the best estimates may be obtained from the unweighted 
estimates with a correction. The structure of the analysis of variance 
and estimation of effects and interactions may be illustrated by the 
following analysis table (Table 4) which can be constructed. For each 


TABLE 4 ANALYSIS TABLE 
EFFECTS AND INTERACTIONS IN EACH REPLICATE 


Replicate Replicate Replicate Weighted 
x Y 4 Mean 


to 


Q 

| 
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replicate the totals of plots for which the level of pseudo-factor a is 
zero, for which it is one, and for which it is two, may be evaluated and 
denoted by A’o, A’;, A’, respectively. Similarly, obtain for b and c the 
totals B’,, B’2, and C’s, and for each interaction, quantities 
like ABC’,, ABC’,, ABC’,. The results may be arranged as in Table 4. 

Asterisks are inserted to denote the effects and interactions which 
are confounded in each replicate. Final estimates of these effects and 
interactions are obtained by forming weighted means of estimates in 
each replicate across the table, using weights w when the effect is un- 
confounded and weights w’ when it is confounded, i.e., when there is an 
asterisk. This table is also useful in the evaluation of the sum of 
squares for blocks adjusted for varieties. Component (a) is obtained 
by taking, for each effect which is confounded in two replicates, the 
differences between the values in the two replicates and obtaining the 
sum of squares between the three differences with the appropriate nu- 
merical divisor. Component (b) is obtained by constructing for each 
of these effects twice the unconfounded values minus the sum of the 
confounded values and evaluating the sum of squares for the differences 
of these three quantities. Component (c) is obtained by evaluating for 
each effect confounded in, one replicate, twice the confounded effect 
minus the sum of the unconfounded effects in the two replicates, and 
evaluating the sum of squares between the three quantities, again yield- 
ing two degrees of freedom for each effect. 

Having obtained estimates of effects and interactions which utilize 
the inter-block information, the adjusted variety totals may be obtained 
by the use of the formula given above. This is computationally a 
tedious process and the value of Yates’ scheme of computation lies in 
the fact that this is simplified by the representation of a cube in two 
dimensions, with the necessity only of adjusting tables involving two 
factors. The method is described here not necessarily as the best com- 
putationally, but as the one which gives clearly the structure of the 
analysis and which may be extended to deal with any number of repli- 
cates for prime power lattices. The great bulk of the work of the 
analysis of variance and calculation of variety means may be dealt 
with simply by the use of punched-card machines. The evaluation of 
the variance of differences between variety means and the discussion of 
the efficiency of the designs are presented in a later section. 

(e) p* system or 4-dimensional lattice in blocks of p. For the pur- 
poses of illustration let the pseudo-factors be denoted by a, b, c, d, each 
having 3 levels. As indicated in the construction of Table 1 for the 3° 
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system, it is necessary in the enumeration of possible schemes of con- 
founding only to designate effects or interactions which may be used to 
generate the whole group of effects and interactions which are con- 


founded. The 40 (equals p*+p*?+p+1) possible schemes are desig- 
nated in Table 5. 


TABLE 5 


ScHEMES OF CONFOUNDING FOR 34 SYSTEM 
IN BLocks or 3 PLots 


Generators of group of confounded interactions 


Number Generators Number Generators 
1 A, B,C 21. B, AC?, AD 
2. A, B, D 22. B, AC?, AD? 
3. A, B, CD 23. AB, C,D 
4. A, B, CD 24. AB, C, AD 
5. 4, C, D 25. AB, C, AD? 
6. A, C, BD 26.._ AB, AC, D 
z. A, C, BD? 27. AB, AC, AD 
8. A, BC, D 28 AB, AC, AD? 
9. A, BC, BD 29 AB, AC, D 
10. A, BC, BD? 30. AB, AC?, AD 
11. A, BC?, D 31 AB, AC?, AD? 
12. A, BC, BD 32. AB?, C, D 
13. A, BC?, 33. AB?, AD 
14. B,C, D 34 AB?, C, AD? 
15. B,C, AD 35 AB?, AC, D 
16. B, C, AD? 36. AB?, AC, AD 
17. B, AC, D 37 AB?, AC, AD? 
18. B, AC, AD 38 AB?2, AC?, D 
19. B, AC, AD? 39 AB?, AC?, AD 
20. B, AC2, D 40. AB?, AC2, AD 


The minimum number of replicates which may be used is four, X, 
Y, Z, W, say, the confounding being 


replicate X: A, B, C and their interactions 


y: A, B, D 
Z: A, C, D 
W: B, C, D “ce 


Each main effect is then confounded in three of the four replicates. 
The structure of the analysis of variance is given in Table 6. 

Blocks component (a) will consist of the comparison between replicates 

of effects which are confounded in those replicates. Thus each of the 
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effects A, B, C, D yields 4 degrees of freedom, with a total of 16. Those 
for A are the interaction of the table 


(A)o (A); (A)s 
Z: - 


TABLE 6 


STRUCTURE OF ANALYSIS OF VARIANCE FOR p* VARIETIES IN 
BLocgs of p PLoTs WITH 4 REPLICATES 


DP Expectation of Mean 
Square 
General p=3 
Replicates 3 3 
Blocks Component (a) 8(p-1) 16 os +p on 
(b) 6(p—1)3 24 os +p ov? 
(d) 6(p—1)2 24 +5? 
(e) 4(p-1)? 82: 
Varieties ignoring blocks pel 80 
Error (intra-block) 8p4p*+1 136 oe 
Total 3238 


Blocks component (b) arises from interactions which are confounded 
in 2 of the 4 replicates. Each of the interactions AB, AB*, AC, AC?, 
AD, AD’, BC, BC?, BD, BD*, CD, CD? yield 2 degrees of freedom with 
a total of 24 degrees of freedom. That for AB is obtained as the inter- 
action of the table. 


(AB), (AB), (AB); 
replicate X: - 


Y: - 


Component (c) arises from the comparison of the mean confounded 
effects and mean unconfounded effects and has 8 degrees of freedom in 
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all: for A, for example, the following table may be constructed. 


Ay A, A, 

replicate X: - 

¥: - 

Z: - 

W: - 


The comparison of the three quantities in the last row of the table yields 
2 degrees of freedom. Component (d) arises in the same way from the 
2-factor interactions and has in all 24 (=12x2) degrees of freedom. 
Finally component (e) arises from the interactions which are con- 
founded in only one replicate. For each, the effect in the replicate in 
which it is confounded is considered and also the effect in the other 
three replicates. From these a comparison yielding 2 degrees of free- 
dom is obtained. There are 16 such interactions yielding 32 degrees of 
freedom in all. The expectation of the mean squares in terms of o»’ 
and o;,? as previously defined is given in the last column of the table. 
The sum of squares for the several block components may be pooled and 
used with the mean square for error in order to estimate o,? and o;? and 
hence w nm the weight of intra-block estimates, and w’ -—_~——, : 5, the 
+ Pon 
weight of inter-block estimates. 

From this stage the analysis could be made in the same way as 
described for 3° system and there is no need to describe it here. Not 
more than 4 replicates are dealt with here, as this extension presents no 
difficulty. 

(f) p® system or 5-dimensional lattice in blocks of p. This case is 
considered as a further illustration of the principles, again with p = 3. 
Let the factors be a, b, c, d, e. The minimum number of replicates is 
5 and the confounding would then consist of the following 


replicate V: A, B, C, D and all their interactions 


XxX: A, B, D, E “cc 
y: A, C, D, E 


The portion of the analysis of variance dealing with blocks with 
400(=5 x 80) degrees of freedom is given in Table 7. 
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TABLE 7 


PARTITION OF BLOCK DEGREES OF FREEDOM FOR p5 VARIETIES 
IN BLOCKS OF p WITH 5 REPLICATES 


DF Expectation of Mean 
General case . p=3 Square 
Comparison between replicates of 
effects or interactions which are 
confounded in these replicates. 
Main effects 15(p-1) 30 
Two-factor interactions 20(p—1)2 80 + 
Three-factor interactions 10(p-1)* 80 
Comparisons between mean of 
confounded effects and mean of | 
unconfounded effects. 
2 
Main effects 5(p-1) 10 5 pov® 
2 
Two-factor interactions 10(p—1)* 40 Por? 
3 
Three-factor interactions 10(p-1)* 80 out 
Four-factor interactions 5(p-1)* 80 + 


The entries in the last column in the second part of the above table 
are obtained by use of the general rule that the expectation of the mean 
square derived from the comparison, when the effect or interaction is 
confounded in n, replicates and unconfounded in n, replicates, where 
blocks of p plots are used, is 


Ny + 


V. MEAN VARIANCE OF ADJUSTED VARIETAL COMPARISONS 


Any comparison between a pair of varieties may be expressed in 
terms of the effects and interactions, using the notation above. Thus in 
the 3° system for example, the comparison of varieties 000 and 001 may 
be expressed as follows: 


Yooo Yoo. = [(C)o— (C)1] + [(AC)o— (AC),] + [ (AC?) (AC?) 2] 
+ [(BC).- (BC),] + [(BC*).— (BC?)2] 
+ [(ABC),-(ABC),] + [(ABC?),— 
+ [(AB?C),—- (AB?C),] + [(AB?C?)5— (AB?C?) 2] 


Furthermore, the contrasts in the square brackets are independent of 
each other, so that their variance may be added to give the variance of 
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the variety difference. If three replicates have been used confounding 
A, B, AB, AB? in the first, A, C, AC, AC? in the second, and B, C, BC, 
BC? in the third replicate, then the variances of the effects and inter- 
actions comparisons are as follows: 


2 
effects A, B, C have a variance + interactions AB, AB?, AC, 


AC*, BC, BC* have a variance of ——";—~, and interactions ABC, 
9(w + 2w’) 


ABC?, AB*C, AB*C* have a variance of —-—~. 
9(3w) 


The variance of the above variety comparison is then 


2 1 4 4 
+ * «| ; 

The possible 26 comparisons of a variety with a chosen one will have 
differing accuracies depending upon whether the comparison involves 
one, two, or three of the main effects, but for practical purposes a single 
measure of accuracy for all comparisons, namely the mean variance of 
all the independent comparisons, is sufficient. This variance may be 
obtained very easily for prime lattice designs with any number of 
replicates by noting that each effect or interaction enters into the same 
number of the comparisons of one variety with the remainder. In the 
particular case of the 3° system, each of the 26 independent comparisons 
involve nine effects or interactions : 4 effects or interactions do not enter 
into any particular comparison—in the above case, these are A, B, AB 
and AB*. Each effect or interaction comparison enters into 18 of the 
total of 26 comparisons, so that the mean variance for all comparisons is 
equal to 


[V(A) +V(B) + V(C) + V(AB) + V(AB*) + V(AC) + V(AC?) 


_+V(BC) + V(BC*) + V(ABC) + V(ABC?) + V(AB*C) 
+ V(AB?C?)]. 
where V(A) for example is the variance of each of Ao, Ai, Az appro- 
priate for comparisons with each other. The mean variance is there- 
fore equal to 


9 ° 26|w+2w Qwiw 3w 


= 
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The formula may be generalized immediately to the case of a p" system 
or the n-dimensional lattice in blocks of p plots with r replicates; the 
mean variance is 


where 
m, is the number of effects or interactions confounded in (r—1) replicates, 


ete. 
and n, is the number of effects or interactions not confounded in any 
replicate. The ; must satisfy the two relations, 


(the total number of effects and 
p-1 interactions), 
and 
r(p*-?-1) 
—-, (the total number 
p-1 of degrees of 

freedom for 
blocks within 
replicates di- 
vided by (p 
1)) 


Thus with the 3° system with 4 replicates (p = 3, n=3, r= 4) and con- 
founding as follows: 


replicate 1: A, B, AB, AB? 
” 2:A, C, AC, AC? 
sa 3:B, C, BC, BC? 
= 4: AB, AC, BC?, and AB?C?, 


it is seen that 1n,=0, 
m,=6, namely A, B, C, AB, AC, BC? 
m,=4, namely AB?, AC?, BC, 
m,=3, namely ABC, AB?C, ABC?, 
and the mean variance is 


13, 2w+2w’ 3wiw’ 4w]- 


A frequently used design will be the p* design in blocks of p plots 
with » replicates confounding all main effects but one and all their 


interactions in each replicate. This may be called the n-dimensional 
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lattice for p" varieties with n replicates in blocks of p plots. The mean 
variance in this case is 


n 
p"-1 |w+(n-1)w’ Qw+(n-2)w’ 


(p—1)?nC, (p—1)*-*] 
3w+(n—3)w | 


n! 
where ,C, equals si(n—s)! 


VI. EFFICIENCY OF DESIGNS 


In studying the efficiency of lattice designs, 3 comparisons are of 
particular interest : 

(a) comparison with randomized blocks which make no use of con- 

founding, 

(b) comparison of information per plot with varying numbers 

and types of replicates, 

(c) comparison of designs with varying sizes of block and hence 

varying amounts of confounding. 

The first type of comparison is in general of only academic interest, 
as it is known that with the recovery of inter-block information the 
lattice design is at least as efficient as the completely randomized design, 
apart from a trivial loss of information due to inaccuracies in weighting 
the confounded and unconfounded effects and interactions. If the- 
costs of analysis of the two types of experiment differed by an amount 
appreciable in comparison with the field and other work of the experi- 
ments (which are approximately the same in both cases), it would be 
necessary to make this comparison. 

The other types of comparison are of importance in the design of 
particular experiments. Comparison type (a) above may be regarded 
as a special case of comparison type (c), but for present purposes, it is 
considered separately, as comparison (c) will be discussed in a succeed- 
ing paper dealing with blocks of size which is a power of p. 

In considering the efficiency of designs, no account is taken of the 
loss in information which results from inaccuracies in the weights. 
The examination of this point by Yates [8,9] and Cochran [2] indi- 
cates that it is of trivial importance. 

(a) Comparison with randomized blocks. The mean variance per 
comparison in the general case of a p" design in blocks of p plots with 
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r replicates has been given above. If the experiment had been ar- 
ranged so that each replicate consisted of a randomized block of p" 
plots, the average error variance would be 


p*—p*-} 
w 


as (p"-!—1) degrees of freedom would have a variance of =, and the 


remaining (p*-?—1) = p"— p"-' degrees of freedom would have 
a variance of 1/w. The mean variance per comparison would be this 
quantity multiplied by 2/r. 

The efficiency of the lattice design relative to randomized complete 
blocks is then 


1 
1 w w 


r(p—1) 


w+ (r—1)w’ (r—2)w’* 


If, for example, w/w’ equals unity, that is, inter-block and intra-block 
information are of equal accuracy, the above expression reduces to 


r(p-1) 
(m+ M+... 
r 
which, as 
... 


is equal to unity. The relative efficiency may be easily computed from 
the above formula, given a system of confounding and a value for w/w’. 

Table 8 deals with only the particular case of the p" lattice in 
blocks of p plots with the least possible number of replicates, namely n. 
In this case, the relative efficiency is 


(9 
l @ -1) 


2w/w’ + (n-2) nw/w’ | 


Most frequent use of lattice designs is made by plant breeders who 
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TABLE 8 


RELATIVE EFFICIENCY oF n-DIMENSIONAL Lattice Designs For p* VARIETIES 


COMPARED TO RANDOMIZED COMPLETE BLOCKS (” REPLICATES 
WITH INCOMPLETE BLOCK SIZE=p) 


0/w’ 

p " p* a 2 3 4 5 6 7 8 9 10 
2 2 4 100 109 125 148 162 181 200 220 239 259 
2 8 8 100 110 127 146 165 185 205 225 245 266 
2 4 16 100 "110 127 146 165 185 205 225 245 266 
2 5 32 100 110 128 147 166 186 206 226 246 266 
2 6 64 100 111 128 147 167 #187 208 228 248 269 
2 7 128 100 111 128 148 168 189 208 230 251 271 
2 8 256 100 111 129 149 170 190 211 232 253 274 
2 9 512 100 111 130 150 170 192 213 234 256 277 
2 10 1024 100 111 130 150 171 193 214 236 258 280 
2 11 2048 100 111 130 151 172 194 216 238 260 282 
8 2 9 100 107 120 135 150 166 182 198 214 231 
3 8 27 100 108 128 189 156 178 191 208 226 244 
3 4 81 100 109 124 141 159 177 195 213 232 250 
8 5 243 100 109 125 142 161 179 198 217 236 255 
8 6 729 100 110 126 144 162 182 201 220 240 259 
8 7 2187 100 110 126 145 164 183 203 228 242 262 
5 2 25 100 105 114 125 186 148 160 172 184 196 
5 3 125 100 106 117 129 #142 155 168 182 196 209 
5 4 625 100 106 118 181 144 158 172 186 200 215 
5 5 8125 100 107 119 132 146 160 174 189 204 219 
7 2 49 100 104 111 120 129 188 147 157 #167 176 
7 8 8438 100 105 113 123 183 143 154 165 176 186 
7 4 2401 100 105 114 124 135 146 157 168 179 190 

11 2 121 “100 108 108 114 120 127 183 140 147 154 

11 8 1381 100 108 109 116 128 180 138 145 153 160 

18 2 169 100 103 107 112 1f7 123 129 185 141 147 

18 8 2197 100 108 108 114 120 126 188 139 146 152 


often have no particular number of varieties to compare. For ex- 
ample, the number of varieties may be of the order of 500 or 600; in 
order to use an efficient experimental design, the geneticist will con- 
strict this number, and if he is considering prime-power lattice 
designs he may use 2°(=512), or 5‘(=625), or 23*(=529) varieties. 
To determine which choice he should make depends on several consider- 
ations among which is the relative efficiency of the two designs. A 
clear-cut answer is not possible because the relative efficiency depends 
on the ratio w/w’ in the two cases. 

(b) Comparison of information per plot with different numbers of 
replicates. Earlier in the paper consideration was given to the use of 
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any number of replicates equal to or greater than the number of pseudo- 
factors (which have been denoted by ”). Particular attention was 
given to the 3° system when it was shown that the 13 possible schemes 
of confounding could be arranged, intuitively at least, in an order such 
that the first r schemes would be the best to use if r replicates were to 
be used. Here this problem is considered in detail with reference to 
the 3° and 3* systems. 

For each number of replicates the efficiency of the design relative to 
complete randomized blocks, may be evaluated, and for particular 
values of w/w’ the relative efficiency for varying numbers of replicates 
may be compared. This comparison is equivalent to comparing the 
information per plot of the various designs, as the information per plot 
with complete randomized blocks is independent of the number of 
replicates. 

For the 3 replicates given by the first 3 rows of Table 2, the relative 
efficiency is 


[18/( w/w’) +8] 


6 4 | 
w/w’ +2 2w/w’ +1 3w/w’ 
For the 4 replicates given by the first 4 rows of Table 2, the relative 
efficiency is 
[18/(w/w’) +8] 


| 
2w/w’'+2 3w/w'+1 4w/w’ 
and so on. 


The values of the relative efficiency for numbers of replicates from 
3 to 13 as specified by Table 2 are given in Table 9. 

In the literature on the 3° system, it is reeommended that the num- 
ber of replicates be a multiple of 3, consisting of the first three repli- 
cates repeated the requisite number of times. From Table 9 a com- 
parison may be made of the efficiency of an experiment of this type with 
6 replicates using different systems of confounding with that of an 
experiment in which the confounding is the same in pairs of replicates. 
Of course, if w/w’ is equal to unity, nothing is to be gained from the 
use of different replicates but for moderate values of w/w’ the efficiency 


is around 10 per cent greater for the design with different systems of 
confounding. 
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TABLE 9 


PERCENTAGE EFFICIENCIES OF THE 33 LATTICE DESIGN FOR VARIOUS 
NuMBERS OF REPLICATES TO RANDOMIZED COMPLETE BLOCKS , 


Systems from Table 1 Values of w/w’ 

1 2 3 4 5 6 7 8 9 10 
1,2,5 100 108 123 189 156 173 191 208 226 244 
1,2,5,9 100 108 125 143 161 180 199 218 238 257 
1,2,5,9,13 100 110 126 144 168 182 202 222 241 261 
1,2,5,9,13,12 100 110 127 145 164 184 204 224 244 264 
1,2,5,9,13,12,10 100 110 127 146 165 185 205 226 246 266 
1,2,5,9,13,12,10,11 100 110 127 146 166 186 206 227 247 268 
1,2,5,9,18,12,10,11,8 100 110 128 147 167 187 207 228 248 269 
1,2,5,9,13,12,10,11,8,3 100 110 128 147 167 188 208 229 250 270 
1,2,5,9,13,12.10,11,8,3,6 100 111 128 148 168 188 209 229 250 271 
1,2,5,9,13,12,10,11,8,3,6,7 100 111 128 148 168 188 208 230 251 272 
1,2,5,9,13,12,10,11,8,3,6,7,4 100 111 128 148 168 189 210 230 252 273 


For the 3‘ system, it was not considered worth while to make as com- 
plete an examination and in Table 10, the relative efficiency of 4, 5, 6 


and 40 replicates are compared. 


TABLE 10 


RELATIVE INFORMATION PER PLOT WITH 4, 5, 6, AND 
40 REPLICATES OF 34 SYSTEM 


w/w’ 
Sets 
1 2 3 4 ° 5 6 7 8 9 10 
1,2,5,14 100 109 124 141 159 177 195 213 232 250 
1,2,5,14,27 100 108 8125 143 162 181 200 220 289 #«=259 
1,2,5,14,27,40 100 110 126 144 163 183 202 222 242 262 
All 40 sets 100 111 129 149 170 191 2138 234 256 278 


Again, when w/w’ is small, there is little gain from using more than 
4 replicates (or multiples of 4 replicates). When w/w’ is large, gains 
of up to 10 per cent are obtained. This, of course, supposes that the 
same size and shape of plots, blocks, and replicates are used in the 


possible cases. 
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QUERIES 


QUERY. In our investigations of processes of freezing vegetables, we 
56 prepared samples from the same lot of raw material and subjected 
_ half to process A and half to a second process B. In order to learn 
if the eventual eating quality is different, and if so which product is 
preferred, the samples were submitted to graders on two days: one 
day AAB was submitted, another, ABB. The order of presentation 
was randomized each day. The grades were asked if they could detect 
any difference between the samples; and if so, which were identical. 
They were also to indicate in each trio the sample or samples they liked 
best. 

Among 150 graders who were present both days, 17 correctly iden- 
tified duplicates. Of these 17, one preferred A both days, 7 pre- 
ferred B both days, and 8 preferred A one day and B the other. One 
had no preference either day. 

Now we can state one of our questions. What is the value of r in 


the formula, 


(a—rb)? 
r(a+b) 


when analyzing the data for those persons who identified duplicates 
both days? We think the grader had one chance in 9 to identify dupli- 
eates on both days, therefore, r = 8. 

Now we will approach our second question which relates to the 
preference data for both days. These data obviously support the con- 
clusion that there is no important difference between the samples as 
far as these graders are concerned. But what if our data had been 
slightly different, say that 30 people had identified duplicates both 
days but of these 30, 15 preferred A one day but B the other day. 
Would we be justified in concluding from the 30 that the samples were 
different? On the other hand, maybe we should think of the difference 
data and the preference data as two separate categories and analyze 
each separately. What do you think? 


Under the hypothesis that the graders can detect no differ- 
ANSWEE. ence between the products, the probability of correct as- 
sociation each day is 1/3 so that the probability of success 
both days is 1/9; you are right in taking r= 8 for substitution in the 
formula quoted provided a=133 and b=17. As you indicate, the 
evidence against the hypothesis is trivial. 
80 
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The preference data suggest a modification in this conclusion. 
Clearly the person who had no preference either day provides no in- 
formation on this point. Of the other 16, half had the same preference 
on both days while the other half changed. So far the results are 
exactly in accord with the hypothesis that A and B are the same. On 
the other hand, of the 8 who had consistent preferences, 7 chose B. 
This is an indication, though the evidence is not strong, that there may 
be a margin of preference in favor of B in the population of graders 
who have consistent preferences. 

I am not sure that I understand the meaning of your second ques- 
tion. Certainly, if it is among 150 graders that there are 30 successful 
identifications of duplicates on successive days, then there is strong evi- 
dence of a difference between A and B. As to preference, one may 
still set up the hypothesis that this is entirely random. If it is, the 
following statements of preference are equally likely : 


Day 1 Day2 
A A 
B B 
A B 
B A 


I suppose you mean that 15 of the 30 graders stated their preferences 
as either AB or BA. This is the number to be expected under the 
hypothesis set up. To this point, the preference results are similar to 
those in your first question. But no information is given about the 
consistent preferences. If they are somewhat evenly divided, the 
hypothesis of random choice need not be rejected. But if most of 
_ them preferred B both days, the results of your first experiment are 
accentuated. 

As for the third question, it seems to me that each part of the 
experiment contributes information on the other part. It may be that 
progress can be made toward segregating the insensitive graders from 
the sensitive, especially if it is found that these qualities are general 
rather than specific. If the generally insensitive are gradually elimi- 
nated, there is a saving of time and money. Both parts of the experi- 
ment are useful in making the separation. Those who both identify 
duplicates and have consistent preferences are promising candidates 
for further testing. 

Grorcs W. SNEDECOR 
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THE BIOMETRIC SOCIETY 


The Biometric Society, a very flourishing infant, is gaining daily. 
Any count of members will be far out of date by the time of publica- 
tion. As of March 1, although the majority of the 400 members 
were from the Continental United States, 35 from other countries 
represented beginnings in Argentina, Australia, Canada, Czechoslo- 
vakia, Denmark, England, France, Hawaii, India, Italy, the Nether- 
lands, Norway, Sweden and Switzerland. 

More interesting than the number of members are the fascinating 
‘mathematical and statistical aspects of biology’’ evidenced by mem- 
bers interested and working in everything from aviation medicine to 
the several facets of zoology. As if the millenium had arrived and 
lions were in fact lying down with lambs, animal and plant breeding 
geneticists, serologists, nutritionists, psychiatrists, ornithologists, phar- 
macologists, biological assayists, sampling and quality control statis- 
ticians, pomologists, herpetologists, cancer and radium biologists, for- 
esters and geologists seem to be brothers and sisters under the skin 
with the biophysicists, mathematical and vital statisticians, at least in 
the cosmopolitan Biometric Society. 

The annual business meting for formation of the Rain North 
American Region was held in two parts, the first at the Stevens Hotel 
in Chicago, December 27, and the second at the Commodore Hotel in 
New York, December 30, 1947. 

Geoffrey Beall presided at the Chicago meeting which took the fol- 
lowing action: 

1. After consideration the proposed By-laws for the Region 
were adopted with a few minor changes. 

2. A slate of officers for 1948, was sent to the Nominating 
Committee in New York to be considered with other nominations. 

3. A Resolution of Appreciation to the American Statistical 
Association was indorsed unanimously : 


The Eastern North American Region of the Biometric So- 
ciety, at its first business meeting, wishes to record the indebt- 
edness of American biometricians to the American Statistical 
Association. The continuing support of the Association has 
made it possible for biometrics to develop into a recognized 
field in America. All biometricians owe much to the Associa- 
tion for its organization of the Biometrics Section, and the 
82 
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founding of Biometrics. It is the hope of the Eastern North 
American Region to establish and maintain the most cordial 
relations with the American Statistical Association. 


4. A motion authorizing the Regional Committee to affiliate 
as soon as feasible with the ASA, AAAS and NRC Division on 
Biology and Agriculture carried unanimously. 

5. A motion that meetings in 1948 be arranged jointly by the 
Biometrics Section of the ASA and the Eastern North American 
Region of the Society was carried unanimously. 

After some discussion of desirable meetings, it was agreed that 
everyone present approved such cooperative endeavors but the 
details should necessarily be left to the Regional Committee. 

6. After considerable debate on the manner in which the spirit 
might be best implemented, it was unanimously agreed that ‘‘the 
meeting in Chicago gives proxy to the acting secretary of the New 
York meeting to make any minor modifications he may see fit in 
action taken in Chicago and to make decisions subject to confirma- 
tion or reversal at the next annual business meeting.’’ 

C. I. Bliss, acting for D. B. DeLury, presided at the New York 
meeting, which by virtue of a common agenda had the objective 
of completing action started at Chicago. 

1. After discussion of the modified, proposed by-laws, and a few 
voted changes, a motion was carried to accept the following: 

As a subdivision of the Biometric Society, the Eastern North 
American Region is governed by the Constitution of the Society and 
the following By-laws. 

1) Scope. The Region shall be concerned with biometrical 
activities, ag defined in Article 1 of the Constitution, in that part 
of the United States and Canada lying east of approximately 104° 
longitude. 

2) Affikations. By action of the Regional Committee, ap- 
proved by the Council of the Society, the Region may affiliate 
itself with national organizations in the fields of biology, mathe- 
matics, and statistics. Such affiliation shall be reported to the 
membership at the next annual meeting for confirmation. 

3) Regional Committee. The Vice-President of the Society 
for the Region, who shall serve as its chairman, the Regional Sec- 
retary-Treasurer and 6 ordinary members shall constitute the 
Regional Committee, all to be elected by the Council of the Bio- 
metric Society. The Regional Committee shall be responsible 
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for the affairs of the Region between annual meetings. Each an- 
nual meeting shall submit nominations for a Vice-President, a Sec- 
retary-Treasurer, and two ordinary members for 3 year terms, to 
the Council. No ordinary member shall be nominated to succeed 
himself. 

4) Regional Advisory Board. There shall be a Regional Ad- 
visory Board, composed of the Regional Committee, members 
designated to represent the Region in other organizations, chair- 
men and other designated representatives of committees, and such 
other members as the Regional Committee shall appoint. The 
Regional Advisory Board shall advise the Regional Committee on 
matters of policy and the appointment of committees, and shall, 
through its members, help to keep the Regional Committee in touch 
with the interests and wishes of all members of the Region. 

5) Annual Meeting. There shall be an annual meeting of the 
Region at a time and place determined by the Regional Committee. 

6) Associate Members. Members of organizations with which 
the Region is formally or informally affiliated who subscribe to 
Biometrics shall be associate members of the Region. They shall 
have all privileges except those of voting and hoiding office. 

7) Program Commésttees. The Vice-President shall appoint 
standing program committees with rotating membership to arrange 
meetings in connection with those of specific organizations or 
groups of organizations. Each such committee shall report its 
tentative plans for the coming year at the annual meeting. Each 
independent meeting of the Region shall be arranged by a special 
program committee appointed for thé purpose. 

8) Nominating Committee. At the annual meeting, the Vice- 
President shall appoint a Nominating Committee of four mem- 
bers, which shall present its nominees for the Regional Committee 
to the next annual meeting. 

9) Amendment. These By-laws may be amended by a-two- 
thirds vote of those attending the annual meeting. Each pro- 
posed amendment shall be placed on the agenda. 

2. The Resolution of Appreciation to the American Statistical 
Association, previously indorsed in Chicago, was indorsed unani- 
mously. 

3. A motion was carried authorizing the Regional Committee to | 
affiliate as soon as’feasible with the ASA, the AAAS and NRC Division 
on Biology and Agriculture. . 
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4. The Chicago motion that meetings in 1948 be arranged jointly 
by Section and Region was confirmed. 

5. The Nominating Committee (H. Muench, L. Knudsen and C. 
Eisenhart) presented a slate of officers. After considerable discussion 
and nominations from the floor, the following officers for 1948 were 
elected for Council confirmation : 


Vice-President : Charles P. Winsor 
Secretary-Treasurer: John H. Watkins 
Regional Committee: 


1948-1950: N. Rashevsky ; H. D. Landahl 
1948-1949 : Margaret Merrell; Phillip J. Rulon 
1948 : E. J. deBeer; A. E. Brandt 


The British Region, which has twenty charter members, held an 
organization meeting in London on January 28. A provisional com- 
mittee (J. W. Trevan, E. C. Fieller, R. A. Fisher, J. B. S. Haldane, 
K. Mather, Eric Wood) was appointed to draft regional rules for 
submission to a later inaugural meeting, and to draw up a list of those 
who should be invited to join. Announcements concerning the Society 
are being sent to Nature and other suitable journals. 

E. A. Cornish, organizing the Australian Region, has sent relevant 
material to everyone in Australia likely to be interested and does not 
think there will be any difficulty at all in setting up the region. Helen 
Turner is acting as Secretary-Treasurer. Cornish has forwarded a list 
of twenty-one who wish to become charter members in the Society. 

A. Buzzati-Traverso has forwarded to the Seeretary’s office a list 
of nine colleagues who have asked to become members of the Society. 
He is getting in touch with Georges Teissier with a proposal that the 
Italian members join the French Region. 
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ANNUAL MEETING OF THE BIOMETRICS 
SECTION 


The annual business meeting of the Biometrics Section was held at 
the Commodore Hotel, New York City, at 2 P.M., December 30, 1947. 
C. I. Bliss, presiding in the absence of the Section Chairman, reported 
on the First International Biometric Conference at Woods Hole, Mas- 
sachusetts. The Proceedings of this conference appeared in the De- 
cember issue of Biometrics. 

The secretary, H. W. Norton, reported that the Section had grown 
to 1040 members and 240 associate members. 

The presiding chairman presented a report of the joint business 
meeting with the Biometrics Society, held at the Stevens Hotel in Chi- 
cago on December 27, 1947. This report, pertaining principally to the 
formation of the Eastern North American Region of the Biometric 
Society, is contained in the report of the Society elsewhere in this issue. 

Gertrude Cox, chairman of the Editorial Committee of Biometrics, 
reported-on the increase in the annual subscription from one to two 
dollars and discussed the revised format and increased size of the 
journal. 

A joint meeting with the Pharmacological Society in Atlantic City 
on March 15-19, 1948 was announced, and a motion was carried that 
meetings during the forthcoming year be arranged jointly with the 
Eastern North Ameriean Region of the Biometric Society. 

The Nominating Committee, C. P. Winsor, chairman, A. E. Brandt, 
P. T. Bruyere, Lila F. Knudsen, Hugo Muench and Churchill Eisen- 
hart, presented a slate of officers for 1948, who were elected unani- 
mously. The new officers are: Chairman, Joseph Berkson; Secretary, 
John H. Watkins ; Section Committee, D. J. Finney, Margaret Merrell, 
P. J. Rulon and J. W. Tukey ; Editorial Committee, Gertrude M. Cox, 
chairman, C. I. Bliss, W. G. Cochran, Churchill Eisenhart, J. W. Fer- 
tig, H. C. Fryer, Horace Norton, A. M. Mood, G. W. Snedecor, and 
Jane Worcester. 
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STATISTICAL SUMMER SESSIONS AT THE UNIVERSITY 
OF CALIFORNIA 


Berkeley, California 


Following the encouraging experience of last year the University 
of California offers statistical programs in two Summer Sessions of 
1948. The teaching staff is as follows: 

Raj Chandra Bose, Professor of the University of Calcutta, 

India. 

Miss Evelyn Fiz, Lecturer at the University of California, Ber- 

keley. 

Erich L. Lehmann, Assistant Professor of the University of 

California, Berkeley. 

Michel Loeve, Reader at the University of London, England. 

Jerzy Neyman, Professor of the University of California, Ber- 

keley. 

Abraham Wald, Professor of Columbia University, New York. 


Courses in statistics are offered on both the graduate and the under- 
graduate levels. The graduate courses, all given during the first Sum- 
mer Session, are meant primarily for students who either have already 
obtained their Ph.D. degree or are working towards it. Therefore, 
apart from formal classes, it is proposed to hold extensive seminars in 
which the work of students will be discussed. No specific prerequisites 
to graduate courses will be required. However, to benefit from the 
courses, the students must be generally familiar with the theory of 
statistics. In addition, course 272 and especially 271 will require a 
reasonable knowledge of the theory of functions. 

There will be two undergraduate courses offered, course $12 during 
the first Summer Session, June 21 to July 31, and course S113 during 
the second Summer Session, August 2 to September 11. Both of these 
courses were recently introduced into the curriculum and are pre- 
requisites to more advanced courses in statistics. They are offered 
during the Summer Sessions for the benefit of students, otherwise 
advanced, who plan to attend more advanced courses in statistics dur- 
ing the fall semester. Besides, course S12 is recommended for students 
who do not intend to specialize in statistics but wish to acquire some 
knowledge of this subject as a part of their general education. 

The Statistical Laboratory will be available for students doing re- 
search. 
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FIRST SUMMER SESSION 
$12. Elements of Probability and 


Statistics. Mr. Lehmann 
271. Random Functions. — Mr. Loeve 
272. Sequential Analysis. Mr. Wald 
273. Design of Experiments. Mr. Bose 
rane Seminar in Theory and Sta- 

tistics. Mr. Loeve, Mr. Wald 
2008 Seminar in Design of Experiments Mr. Bose 
$295. Individual Research. Mr. Bose, Mr. Loeve, 


Mr. Neyman, Mr. Wald 


SECOND SUMMER SESSION 


$113. Second Course in Probability 
and Statistics. Miss Fix 


Other details may be had upon request to the Summer Sessions 
Office, University of California, Berkeley 4, California. 
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NEWS AND NOTES 


George W. Snedecor, President of the American Statistical Asso- 
ciation and Research Professor of Statistics at Iowa State College, will 
be Visiting Research Professor of Statistics at Alabama Polytechnic 
Institute during the Spring Quarter, from March 22 to June 4, 1948. 
He will lecture on Statistical Experimental Design and will be avail- 
able for consultations. The newly formed Statistical Laboratory at 
A. P. I. also will offer a course in Survey Sampling during the Spring 
Quarter to be taught by the Director, T. A. Bancroft. Conferences 
in applied statistics for research workers in the lower southeastern 
states are being scheduled during the time of Mr. Snedecor’s visit .. . 
Horace Norton, formerly with the Washington office of the Weather 
Bureau, reported November 3, for duty with the Atomic Energy Com- 
mission. His new work is in connection with accountability for source 
and fissionable materials .. . Carl F. Kossack keeps us guessing 
where to locate him. At last report he is with the Department of 
Mathematics at Purdue University, Indiana . . . Bliss H. Crandall, 
Director of the Statistical Laboratory, Utah State College, Logan, is 
spending six months at the Institute of Statistics, Raleigh, North 
Carolina. He is visiting classes and assisting with consulting and 
analytical work . . . A number of biologists in the New York Metro- 
politan area have diagnosed themselves as being statistically under- 
nourished and have banded together to form a study group designed 
to alleviate this deficiency. Self-medication consisting of ‘‘t’s’’, re- 
gression coefficients, variances, and assorted statistical treatments are 
administered at bi-weekly meetings. If you are interested in joining 
this clinic, either as an observer or as a patient, get in touch with 
Edwin J. de Beer, The Wellcome Research Laboratories, Tuckahoe 7, 
New York .. . J. I. Northam from the University of Michigan and 
Marily Spanglet from the University of California are now at Kansas 
State College helping teach statistics. H.C. Fryer has a statistical 
program operating at the Agricultural Experiment Station and Kansas 
State College, Manhattan . . . F. J. Anscombe, Rothamsted Experi- 
mental Station, has been appointed lecturer in the Faculty of Mathe- 
matics, University of Cambridge . . . The Calcutta Statistical Asso- 
ciation was started towards the end of 1945 with the object of promot- 
ing the cause of statistics in post-war India in all possible ways. Their 
Bulletin states, ‘‘There is a greater need in our country than elsewhere 
of training up competent statisticians and educating the government 
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and the public on the usefulness of reliable statistics and their dispas- 
sionate analysis in solving the diverse problems of our life.’? Such a 
need is‘not confined to India. In fact, the statisticians of India are 
taking an active part in promoting national and interyational statis- 
tical programs. P. C. Mahalanobis is a member of the Statistical Com- 
mission set up by the Economic and Social Council of the United 
Nations. He is Chairman of the Sub-Commission on Statistical 
Sampling . . . R. C. Bose, Calcutta University, is a Visiting Pro- 
fessor at the University of North Carolina where he is teaching a 
course in Multivariate Analysis. In the Spring Quarter he will give 
a course in Design of Experiments. During the Fall he was at Colum- 
bia University, New York . . . D. Mangeron, Director of the Mathe- 
matical Institute of the Polytechnical High School of Jassy, Roumania, 
writes, ‘‘Let us rejoice that this period of destruction is come to an 
end and the universal desire among scientists to see better international 
scientific contacts in the coming years begin to become a reality.’’ 
They are soliciting original articles of Mathematics, ‘‘Fisics,’’ Chem- 
istry, and Technics for publication and are requesting extracts or re- 
prints of publications destinated eventually to a review in their Bul- 
letin . . . Joseph Carmin, Director, Independent Biological Labora- 
tories, Kefar-Malal, P. O. Ramatayim, Palestine, reports, ‘‘ A building 
was erected at Kefar-Malal, an agricultural settlement... It is a 
pleasure to be able to state that the moving and the rearranging is — 
almost all over and that we are ready now for research work; we are 
able to accommodate a dozen workers with working facilities.’’ 
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